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OLIGONUCLEOTIDE TAGS FOR SORTING 
AND IDENTIFICATION 

5 Field of the Invention - 

The invention relates generally to methods for identifying, sorting, and/or 
tracking molecules, especially polynucleotides, with oligonucleotide tags, and more 
particularly, to a method of sorting and analyzing such tagged polynucleotides by 
specific hybridization of the tags to their complements. 

10 

BACKGROUND 

Specific hybridization of oligonucleotides and their analogs is a fundamental 
process that is employed in a wide variety of research, medical, and industrial 
applications, including the identification of disease-related polynucleotides in 

1 5 diagnostic assays, screening for clones of novel target polynucleotides, identification 
of specific polynucleotides in blots of mixtures of polynucleotides, amplification of 
specific target polynucleotides, therapeutic blocking of inappropriately expressed 
genes, DNA sequencing, and the like, e.g. Sambrook et al, Molecular Cloning: A 
Laboratory Manual, 2nd Edition (Cold Spring Harbor Laboratory, New York, 1989); 

20 Keller and Manak, DNA Probes, 2nd Edition (Stockton Press, New York, 1993); 
Milligan et al, J. Med. Chem., 36: 1923-1937 (1993); Drmanac et al, Science, 260: 
1649-1652 (1993); Bains, J. DNA Sequencing and Mapping, 4: 143-150 (1993). 

Specific hybridization has also been proposed as a method of tracking, 

retrieving, and identifying compounds labeled with oligonucleotide tags. For 

25 example, in multiplex DNA sequencing oligonucleotide tags are used to identify 

electrophoretically separated bands on a gel that consist of DNA fragments generated 
in the same sequencing reaction. In this way, DNA fragments from many sequencing 
reactions are separated on the same lane of a gel which is then blotted with separate 

solid phase materials on which the fragment bands from the separate sequencing 

30 reactions are visualized with oligonucleotide probes that specifically hybridize to 
complementary tags, Church et al, Science, 240: 185-188 (1988). Similar uses of 
oligonucleotide tags have also been proposed for identifying explosives, potential 
pollutants, such as crude oil, and currency for prevention and detection of 
counterfeiting, e.g. reviewed by Dollinger, pages 265-274 in Mullis et al, editors, The 

35 Polymerase Chain Reaction (Birkhauser, Boston, 1994). More recently, systems 

employing oligonucleotide tags have also been proposed as a means of manipulating 
and identifying individual molecules in complex combinatorial chemical libraries, for 
example, as an aid to screening such libraries for drug candidates, Brenner and 
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Lerner, Proc. Natl. Acad. Sci.; 89: 5381-5383 (1992); Alper, Science, 264: 1399-1401 
(1994); and Needels et al, Proc. Natl. Acad. Sci., 90: 10700-10704 (1993). 

The successful implementation of such tagging schemes depends in large part 
on the success in achieving specific hybridization between a tag and its 
5 complementary probe. That is, for an oligonucleotide tag to successfully identify a 
substance, the number of false positive and false negative signals must be minimized. 
Unfortunately, such spurious signals are not uncommon because base pairing and base 
stacking free energies vary widely among nucleotides in a duplex or triplex structure. 
For example, a duplex consisting of a repeated sequence of deoxyadenosine (A) and 

10 thymidine (T) bound to its complement may have less stability than an equal-length 
duplex consisting of a repeated sequence of deoxyguanosine (G) and deoxycytidine 
(C) bound to a partially complementary target containing a mismatch. Thus, if a 
desired compound from a large combinatorial chemical library were tagged with the 
former oligonucleotide, a significant possibility would exist that, under hybridization 

1 5 conditions designed to detect perfectly matched AT-rich duplexes, undesired 

compounds labeled with the GC-rich oligonucleotide-even in a mismatched duplex- 
would be detected along with the perfectly matched duplexes consisting of the AT- 
rich tag. In the molecular tagging system proposed by Brenner et al (cited above), the 
related problem of mis-hybridizations of closely related tags was addressed by 

20 employing a so-called "comma-less" code, which ensures that a probe out of register 
(or frame shifted) with respect to its complementary tag would result in a duplex with 
one or more mismatches for each of its five or more three-base words, or "codons." 

Even though reagents, such as tetramethylarnmonium chloride, are available to 
negate base-specific stability differences of oligonucleotide duplexes, the effect of 

25 such reagents is often limited and their presence can be incompatible with, or render 
more difficult, further manipulations of the selected compounds, e.g. amplification by 
polymerase chain reaction (PCR), or the like. 

Such problems have made the simultaneous use of multiple hybridization 
probes in the analysis of multiple or complex genetic loci, e.g. via multiplex PCR, 

30 reverse dot blotting, or the like, very difficult. As a result, direct sequencing of 

certain loci, e.g. HLA genes, has been promoted as a reliable alternative to indirect 
methods employing specific hybridization for the identification of genotypes, e.g. 
Gyllensten et al, Proc. Natl. Acad. Sci., 85: 7652-7656 (1988). 

The ability to sort cloned and identically tagged DNA fragments onto distinct 

35 solid phase supports would facilitate such sequencing, particularly when coupled with 
a non gel-based sequencing methodology simultaneously applicable to many samples 
in parallel. 
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In view of the above, it would be useful if there were available an 
oligonucleotide-based tagging system which provided a large repertoire of tags, but 
which also minimized the occurrence of false positive and false negative signals 
without the need to employ special reagents for altering natural base pairing and base 
5 stacking free energy differences. Such a tagging system would find applications in 
many areas, including construction and use of combinatorial chemical libraries, large- 
scale mapping and sequencing of DNA, genetic identification, medical diagnostics, 
and the like. 

10 Summary of the Invention 

An object of my invention is to provide a molecular tagging system for 
tracking, retrieving, and identifying compounds. 

Another object of my invention is to provide a method for sorting identical 
molecules, or subclasses of molecules, especially polynucleotides, onto surfaces of 
15 solid phase materials by the specific hybridization of oligonucleotide tags and their 
complements. 

A further object of my invention is to provide a method for analyzing gene 
expression patterns in diseased and normal tissues. 

A still further object of my invention is to provide a system for tagging and 
20 sorting many thousands of fragments, especially randomly overlapping fragments, of 
a target polynucleotide for simultaneous analysis and/or sequencing. 

Another object of my invention is to provide a rapid and reliable method for 
sequencing target polynucleotides having a-length in the range of-a few hundred 
basepairs to several tens of thousands of basepairs. 
25 A further object of my invention is to provide a method for reducing the 

number of separate template preparation steps required in large scale sequencing 
projects employing conventional Sanger-based sequencing techniques. 

My invention achieves these and other objects by providing a method and 

materials for tracking, identifying, and/or sorting classes or subpopulations of 

30 molecules by the use of oligonucleotide tags. An important feature of the invention is 
that the oligonucleotide tags are members of a minimally cross-hybridizing set of 
oligonucleotides. The sequences of oligonucleotides of such a set differ from the 
sequences of every other member of the same set by at least two nucleotides. Thus, 
each member of such a set cannot form a duplex (or triplex) with the complement of 
35 any other member with less than two mismatches. Complements of oligonucleotide 
tags of the invention, referred to herein as "tag complements," may comprise natural 
nucleotides or non-natural nucleotide analogs. Preferably, tag complements are 
attached to solid phase supports. Such oligonucleotide tags when used with their 
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corresponding tag complements provide a means of enhancing specificity of 
hybridization for sorting, tracking, or labeling molecules, especially polynucleotides. 

Minimally cross-hybridizing sets of oligonucleotide tags and tag complements 
may be synthesized either combinatorially or individually depending on the size of the 
5 set desired and the degree to which cross-hybridization is sought to be minimized (or 
stated another way, the degree to which specificity is sought to be enhanced). For 
example, a minimally cross-hybridizing set may consist of a set of individually 
synthesized 10-mer sequences that differ from each other by at least 4 nucleotides, 
such set having a maximum size of 332 (when composed of 3 kinds of nucleotides 
1 0 and counted using a computer program such as disclosed in Appendix Ic). 

Alternatively, a rmmmally cross-hybridizing set of oligonucleotide tags may also be 
assembled combinatorially from subunits which themselves are selected from a 
minimally cross-hybridizing set. For example, a set of minimally cross-hybridizing 
12-mers differing from one another by at least three nucleotides may be synthesized 
1 5 by assembling 3 subunits selected from a set of minimally cross-hybridizing 4-mers 
that each differ from one another by three nucleotides. Such an embodiment gives a 
maximally sized set of 9 3 , or 729, 12-mers. The number 9 is number of 
oligonucleotides listed by the computer program of Appendix la, which assumes, as 
with the 10-mers, that only 3 of the 4 different types of nucleotides are used. The set 
20 is described as "maximal" because the computer programs of Appendices Ia-c provide 
the largest set for a given input (e.g. length, composition, difference in number of 
nucleotides between members). Additional rninimally cross-hybridizing sets may be 
formed from subsets of such calculated sets. 

Oligonucleotide tags may be single stranded and be designed for specific 
25 hybridization to single stranded tag complements by duplex formation or for specific 
hybridization to double stranded tag complements by triplex formation. 
Oligonucleotide tags may also be double stranded and be designed for specific 
hybridization to single stranded tag complements by triplex formation. 

When synthesized combinatorially, an olig o nucleotide tag of the invention 

30 preferably consists of a plurality of subunits, each subunit consisting of an 

oligonucleotide of 3 to 9 nucleotides in length wherein each subunit is selected from 
the same miriimally cross-hybridizing set. In such embodiments, the number of 
oligonucleotide tags available depends on the number of subunits per tag and on the 
length of the subunits. The number is generally much less than the number of all 
35 possible sequences the length of the tag, which for a tag n nucleotides long would be 
4n 

In one aspect of my invention, complements of oligonucleotide tags attached 
to a solid phase support are used to sort polynucleotides from a mixture of 
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polynucleotides each containihg a tag. In this embodiment, complements of the 
oligonucleotide tags are synthesized on the surface of a solid phase support, such as a 
microscopic bead or a specific location on an array of synthesis locations on a single 
support, such that populations of identical sequences are produced in specific regions. 

5 That is, the surface of each support, in the case of a bead, or of each region, in the 
case of an array, is derivatized by only one type of complement which has a particular 
sequence. The population of such beads or regions contains a repertoire of 
complements with distinct sequences. As used herein in reference to oligonucleotide 
tags and tag complements, the term "repertoire" means the set of minimally cross- 

1 0 hybridizing set of oligonucleotides that make up the tags in a particular embodiment 
or the corresponding set of tag complements. 

The polynucleotides to be sorted each have an oligonucleotide tag attached, 
such that different polynucleotides have different tags. As explained more fully 
below, this condition is achieved by employing a repertoire of tags substantially 

1 5 greater than the population of polynucleotides and by taking a sufficiently small 

sample of tagged polynucleotides from the full ensemble of tagged polynucleotides. 
After such sampling, when the populations of supports and polynucleotides are mixed 
under conditions which permit specific hybridization of the oligonucleotide tags with 
their respective complements, identical polynucleotides sort onto particular beads or 

20 regions. The sorted populations of polynucleotides can then be manipulated on the 
solid phase support by micro-biochemical techniques. 

Generally, the method of my invention comprises the following steps: (a) 
attaching an oligonucleotide tag from a repertoireof tags to each molecule in a 
population of molecules (i) such that substantially all different molecules or different 

25 subpopulations of molecules in the population have different oligonucleotide tags 
attached and (ii) such that each oligonucleotide tag from the repertoire is selected 
from the same minimally cross-hybridizing set ; and (b) sorting the molecules of the 
population onto one or more solid phase supports by specifically hybridizing the 

oli gonucleotide tags with their res pective complements attached to such su pports. 

30 An important aspect of my invention is the use of the oligonucleotide tags to 

sort polynucleotides for parallel sequence determination. Preferably, such sequencing 
is carried out by the following steps: (a) generating from the target polynucleotide a 
plurality of fragments that cover the target polynucleotide; (b) attaching an 
oligonucleotide tag from a repertoire of tags to each fragment of the plurality (i) such 

35 substantially all different fragments have different oligonucleotide tags attached and 
(ii) such that each oligonucleotide tag from the repertoire is selected from the same 
rninimally cross-hybridizing set; (c) sorting the fragments onto one or more solid 
phase supports by specifically hybridizing the oligonucleotide tags with their 
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respective complements attached to the solid phase supports; (d) determining the 
nucleotide sequence of a portion of each of the fragments of the plurality, preferably 
by a single-base sequencing methodology as described below; and (e) determining the 
nucleotide sequence of the target polynucleotide by collating the sequences of the 
5 fragments. 

Another important aspect of my invention is the determination of a profile, or 
a frequency distribution, of genes being expressed in a given tissue or cell type, 
wherein each such gene is identified by a portion of its sequence. Preferably, such 
frequency distribution is determined by the following steps: (a) forming a cDNA 
1 0 library from a population of mRNA molecules, each cDNA molecule in the cDNA library 
having an oligonucleotide tag attached, (i) such that substantially all different cDNA 
molecules have different oligonucleotide tags attached and (ii) such that each 
oligonucleotide tag from the repertoire is selected from the same minimally cross- 
hybridizing set; (b) sorting the cDNA molecules by specifically hybridizing the 

1 5 oligonucleotide tags with their respective complements attached to one or more solid phase 
supports; (c) determining the nucleotide sequence of a portion of each of the sorted cDNA 
molecules; and (d) forming a frequency distribution of mRNA molecules from the 
nucleotide sequences of the portions of sorted cDNA molecules. 

My invention overcomes a key deficiency of current methods of tagging or 

20 labeling molecules with oligonucleotides: By coding the sequences of the tags in 
accordance with the invention, the stability of any mismatched duplex or triplex 
between a tag and a complement to another tag is far lower than that of any perfectly 
matched duplex between the tag and its own complement. Thus, the problem of 
incorrect sorting because of mismatch duplexes of GC-rich tags being more stable 

25 than perfectly matched AT-rich tags is eliminated. 

When used in combination with solid phase supports, such as microscopic 
beads, my invention provides a readily automated system for manipulating and sorting 
polynucleotides, particularly useful in large-scale parallel operations, such as large- 
scale DNA sequencing, wherein many target polynucleotides or many segments of a 

30 single target polynucleotide are sequenced and/or analyzed simultaneously. 

Brief Description of the Drawing s 
Figure 1 is a flow chart illustrating a general algorithm for generating 
minimally cross-hybridizing sets. 
35 Figure 2 diagrammatically illustrates an apparatus for carrying out parallel 

operations, such as polynucleotide sequencing, in accordance with the invention. 

Figure 3 illustrates an embodiment for genotyping by sorting ligated probes 
onto a solid phase support. 
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Definitions 

"Complement" or "tag complement" as used herein in reference to 
oligonucleotide tags refers to an oligonucleotide to which a oligonucleotide tag 
specifically hybridizes to form a perfectly matched duplex or triplex. In embodiments 
5 where specific hybridization results in a triplex, the oligonucleotide tag may be 
selected to be either double stranded or single stranded. Thus, where triplexes are 
formed, the term "complement" is meant to encompass either a double stranded 
complement of a single stranded oligonucleotide tag or a single stranded complement 
of a double stranded oligonucleotide tag. 

1 0 The term "oligonucleotide" as used herein includes linear oligomers of natural 

or modified monomers or linkages, including deoxyribonucleosides, ribonucleosides, 
anomeric forms thereof, peptide nucleic acids (PNAs), and the like, capable of 
specifically binding to a target polynucleotide by way of a regular pattern of 
monomer-to-monomer interactions, such as Watson-Crick type of base pairing, base 

1 5 stacking, Hoogsteen or reverse Hoogsteen types of base pairing, or the like: Usually 
monomers are linked by phosphodiester bonds or analogs thereof to form 
oligonucleotides ranging in size from a few monomeric units, e.g. 3-4, to several tens 
of monomeric units. Whenever an oligonucleotide is represented by a sequence of 
letters, such as "ATGCCTG," it will be understood that the nucleotides are in 5'->y 

20 order from left to right and that "A" denotes deoxyadenosine, "C" denotes 

deoxycytidine, "G" denotes deoxyguanosine, and "T" denotes thymidine, unless 
otherwise noted. Analogs of phosphodiester linkages include phosphorothioate, 

phosphorodithioate, phosphoranilidate, phosphoramidate, and the like. Usually 

oligonucleotides of the invention comprise the four natural nucleotides; however, they 

25 may also comprise non-natural nucleotide analogs. It is clear to those skilled in the 
art when oligonucleotides having natural or non-natural nucleotides may be 
employed, e.g. where processing by enzymes is called for, usually oligonucleotides 
consisting of natural nucleotides are required. 
"Perfectly matched" in reference to a duplex means that the poly- or 



30 oligonucleotide strands making up the duplex form a double stranded structure with 
one other such that every nucleotide in each strand undergoes Watson-Crick 
basepairing with a nucleotide in the other strand. The term also comprehends the 
pairing of nucleoside analogs, such as deoxyinosine, nucleosides with 2-aminopurine 
bases, and the like, that may be employed. In reference to a triplex, the term means 

35 that the triplex consists of a perfectly matched duplex and a third strand in which 
every nucleotide undergoes Hoogsteen or reverse Hoogsteen association with a 
basepair of the perfectly matched duplex. Conversely, a "mismatch" in a duplex 
between a tag and an oligonucleotide means that a pair or triplet of nucleotides in the 
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duplex or triplex fails to undergo Watson-Crick and/or Hoogsteen and/or reverse 
Hoogsteen bonding. 

As used herein, "nucleoside" includes the natural nucleosides, including 2'- 
deoxy and 2-hydroxyl forms, e.g. as described in Kornberg and Baker, DNA 
5 Replication, 2nd Ed. (Freeman, San Francisco, 1992). "Analogs" in reference to 
nucleosides includes synthetic nucleosides having modified base moieties and/or 
modified sugar moieties, e.g. described by Scheit, Nucleotide Analogs (John Wiley, 
New York, 1980); Uhlman and Peyman, Chemical Reviews, 90: 543-584 (1990), or 
the like, with the only proviso that they are capable of specific hybridization. Such 

1 0 analogs include synthetic nucleosides designed to enhance binding properties, reduce 
complexity, increase specificity, and the like. 

As used herein "sequence determination" or "determining a nucleotide 
sequence" in reference to polynucleotides includes determination of partial as well as 
full sequence information of the polynucleotide. That is, the term includes sequence 

1 5 comparisons, fingerprinting, and like levels of information about a target 

polynucleotide, as well as the express identification and ordering of nucleosides, 
usually each nucleoside, in a target polynucleotide. The term also includes the 
determination of the identification, ordering, and locations of one, two, or three of the 
four types of nucleotides within a target polynucleotide. For example, in some 

20 embodiments sequence determination may be effected by identifying the ordering and 
locations of a single type of nucleotide, e.g. cytosines, within the target polynucleotide 
"CATCGC ..." so that its sequence is represented as a binary code, e.g. "100101 ... " for 
"C-(not C)-(not C)-C-(not C)-C ... " and the like. 

As used herein, the term "complexity" in reference to a population of 

25 polynucleotides means the number of different species of molecule present in the 
population. 

Detailed Description of the Invention 

The invention provides a method of labeling and sorting molecules, 

30 particularly polynucleotides, by the use of oligonucleotide tags. The oligonucleotide 
tags of the invention belong to minimally cross-hybridizing sets of oligonucleotides. 
Thus, the sequences of any two oligonucleotide tags of a repertoire will never be 
"closer" than differing by two nucleotides. In particular embodiments, sequences of 
any two oligonucleotide tags of a repertoire can be even "further" apart, e.g. by 
35 designing a rninimally cross-hybridizing set such that oligonucleotides cannot form a 
duplex or triplex with the complement of another member of the same set with less 
than three mismatched nucleotides, and so on. In such embodiments, greater 
specificity is achieved, but the total repertoire of tags is smaller. Thus, for tags of a 
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given length, a trade off must be made between the degree of specificity desired and 
the size of repertoire desired. The invention is particularly useful in labeling and 
sorting polynucleotides for parallel operations, such as sequencing, fingerprinting or 
other types of analysis. 

5 

Oligonucleotide Tags and Tag Complements 
The nucleotide sequences of oligonucleotides of a minimally cross-hybridizing 
set are conveniently enumerated by simple computer programs following the general 
algorithm illustrated in Fig. 3, and as exemplified by programs whose source codes 

1 0 are listed in Appendices la and lb. Program minhx of Appendix la computes all 

minimally cross-hybridizing sets having 4-mer subunits composed of three kinds of 
nucleotides. Program tagN of Appendix lb enumerates longer oligonucleotides of a 
minimally cross-hybridizing set. Similar algorithms and computer programs are 
readily written for listing oligonucleotides of minimally cross-hybridizing sets for any 

1 5 embodiment of the invention. Table I below provides guidance as to the size of sets 
of minimally cross-hybridizing oligonucleotides for the indicated lengths and number 
of nucleotide differences. The above computer programs were used to generate the 
numbers. 



20 



Table I 



Oligonucleotid 
e 

Word 
Length 



4 
6 
7 
7 

r- 

8 
8 
9 
10 
10 
11 
18 
18 



Nucleotide 
Difference 
between 
Oligonucleotides 
of Minimally 
Cross- 
Hybridizing Set 



Maximal Size 
of Minimally 

Cross- 
Hybridizing 
Set 



3 
3 
4 
5 

4 

5 
5 
5 
6 
5 
6 
12 



9 
27 
27 

8 

190- 

62 
18 
39 
332 
28 
187 
*25000 
24 



Size of 
Repertoire 
with Four 

Words 



6561 
5.3 x 10 3 
5.3 x 10 3 
4096 

-b3frx-10- 



1.48 x 10' 
1.05 x 10* 
2.31x10° 
1.21xlO IU 
6.15 x 10* 



Size of 
Repertoire with 
Five Words 



5.90 x 10* 
1.43 x 10' 
1.43 x 10' 
3.28 x 10" 
-2.48-x-10 u 
9.16 x 10° 
1.89 x 10° 
9.02 x 10' 

1.72 x 10' 



For some embodiments of the invention, where extremely large repertoires of 
tags are not required, oligonucleotide tags of a minimally cross-hybridizing set may 
25 be separately synthesized. Sets containing several hundred to several thousands, or 
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even several tens of thousands, of oligonucleotides may be synthesized directly by a 
variety of parallel synthesis approaches, e.g. as disclosed in Frank et al, U.S. patent 
4,689,405; Frank et al, Nucleic Acids Research, 11: 4365-4377 (1983); Matson et al, 
Anal. Biochem., 224: 1 10-116 (1995); Fodor et al, International application 
5 PCT/US93/04145; Pease et al, Proc. Natl. Acad. Sci., 9-1 : 5022-5026 (1994); 
Southern et al, J. Biotechnology, 35: 217-227 (1994), Brennan, International 
application PCT/US94/05896; Lashkari et al, Proc. Natl. Acad. Sci., 92: 7912-7915 
(1995); or the like. 

Preferably, oligonucleotide tags of the invention are synthesized 

1 0 combinatorially out of subunits between three and six nucleotides in length and 

selected from the same minimally cross-hybridizing set. For oligonucletides in this 
range, the members of such sets may be enumerated by computer programs based on 
the algorithm of Fig. 3. 

The algorithm of Fig. 3 is implemented by first defining the characteristics of 

1 5 the subunits of the minimally cross-hybridizing set, i.e. length, number of base 

differences between members, and composition, e.g. do they consist of two, three, or 
four kinds of bases. A table M n , n=l, is generated (100) that consists of all possible 
sequences of a given length and composition. An initial subunit S\ is selected and 
compared (120) with successive subunits S[ for i=n+l to the end of the table. 

20 Whenever a successive subunit has the required number of mismatches to be a 

member of the imnimally cross-hybridizing set, it is saved in a new table M n +j (125), 
that also contains subunits previously selected in prior passes through step 120. For 
example, in the first set of comparisons, M2 will contain S] ; in the second set of 
comparisons, M3 will contain S\ and S2; in the third set of comparisons, M4 will 

25 contain S\ 9 S2, and S3; and so on. Similarly, comparisons in table Mj will be 

between Sj and all successive subunits in Mj. Note that each successive table M n +i 
is smaller than its predecessors as subunits are eliminated in successive passes 
through step 130. After every subunit of table M n has been compared (140) the old 

table is replaced b y the new table Mn+i , and the next round of comparisons are 

30 begun. The process stops (1 60) when a table M n is reached that contains no 
successive subunits to compare to the selected subunit S|, i.e. M n =M n +i . 

Preferably, rriinimally cross-hybridizing sets comprise subunits that make 
approximately equivalent contributions to duplex stability as every other subunit in 
the set. In this way, the stability of perfectly matched duplexes between every subunit 

35 and its complement is approximately equal. Guidance for selecting such sets is 

provided by published techniques for selecting optimal PCR primers and calculating 
duplex stabilities, e.g. Rychlik et al, Nucleic Acids Research, 17: 8543-8551 (1989) 
and 18: 6409-6412 (1990); Breslauer et al, Proc. Natl. Acad. Sci., 83: 3746-3750 
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(1986); Wetmur, Crit. Rev. Biochem. Mol. Bio; 26: 227-259 (1991);and the like. 
For shorter tags, e.g. about 30 nucleotides or less, the algorithm described by Rychlik 
and Wetmur is preferred, and for longer tags, e.g. about 30-35 nucleotides or greater, 
an algorithm disclosed by Suggs et al, pages 683-693 in Brown, editor, ICN-UCLA 
5 Symp. Dev. Biol., Vol. 23 (Academic Press, New York; 1981) may be conveniently 
employed. Clearly, the are many approaches available to one skilled in the art for 
designing sets of minimally cross-hybridizing subunits within the scope of the 
invention. For example, to minimize the affects of different base-stacking energies of 
terminal nucleotides when subunits are assembled, subunits may be provided that 

1 0 have the same terminal nucleotides. In this way, when subunits are linked, the sum of 
the base-stacking energies of all the adjoining terminal nucleotides will be the same, 
thereby reducing or eliminating variability in tag melting temperatures. 

A "word" of terminal nucleotides, shown in italic below, may also be added to 
each end of a tag so that a perfect match is always formed between it and a similar 

1 5 terminal 4t word" on any other tag complement. Such an augmented tag would have 
the form: 



w 


w, 


w 2 


••• 


w k _, 


w k 


w 


w 


w,' 


W 2 ' 




W k .,' 


w k - 


w 



where the primed W's indicate complements. With ends of tags always forming 
20 perfectly matched duplexes, all mismatched words will be internal mismatches 

thereby reducing the stability of tag-complement duplexes that otherwise would have 
mismatched words at their ends. It is well known that duplexes with internal 
mismatches are significantly less stable than duplexes with the same mismatch at a 
terminus. 

25 A preferred embodiment of minimally cross-hybridizing sets are those whose 

subunits are made up of three of the four natural nucleotides. As will be discussed 

v_below^tJte-absence_of_one^ 

permits target polynucleotides to be loaded onto solid phase supports by use of the 
5'-+3' exonuclease activity of a DNA polymerase. The following is an exemplary 
30 minimally cross-hybridizing set of subunits each comprising four nucleotides selected 
from the group consisting of A, G, aad T: 



35 
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Table H 

Word: w 1 w 2 w 3 w 4 

Sequence: GATT TGAT TAGA TTTG 



Word: w 5 w 6 w 7 w g 

Sequence: GTAA AGTA ATGT AAAG 

5 

In this set, each member would form a duplex having three mismatched bases with 
the complement of every other member. 

Further exemplary minimally cross-hybridizing sets are listed below in Table 
IIL Clearly, additional sets can be generated by substituting different groups of 
1 0 nucleotides, or by using subsets of known minimally cross-hybridizing sets. 



Table ffl 

Exemplary Minimally Cross-Hybridizing Sets of 4-mer Subunits 

15 



Set 1 


Set 2 


Set 3 


Set 4 


Set 5 


Set 6 


CATT 


ACCC 


AAAC 


AAAG 


AACA 


AACG 


CTAA 


AGGG 


ACCA 


ACCA 


ACAC 


ACAA 


TCAT 


CACG 


AGGG 


AGGC 


AGGG 


AGGC 


ACTA 


CCGA 


CACG 


CACC 


CAAG 


CAAC 


TACA 


CGAC 


CCGC 


CCGG 


CCGC 


CCGG 


TTTC 


GAGC 


CGAA 


CGAA 


CGCA 


CGCA 


ATCT 


GCAG 


GAGA 


GAGA 


GAGA 


GAGA 


AAAC 


GGCA 


GCAG 


GCAC 


GCCG 


GCCC 




AAAA 


GGCC 


GGCG 


GGAC 


GGAG 


Set 7 


Set 8 


Set 9 


Set 10 


Set 11 


Set 12 


AAGA 


AAGC 


AAGG 


ACAG 


ACCG 


ACGA 


ACAC 


ACAA 


ACAA 


AACA 


AAAA 


AAAC 


AGCG 


AGCG 


AGCC 


AGGC 


AGGC 


AGCG 


CAAG 


CAAG 


CAAC 


CAAC 


CACC 


CACA 


CCCA 


CCCC 


CCCG 


CCGA 


CCGA 


CCAG 


CGGC 


CGGA 


CGGA 


CGCG 


CGAG 


CGGC 


GACC 


GACA 


GACA 


GAGG 


GAGG 


GAGG 


GCGG 


GCGG 


GCGC 


GCCC 


GCAC 


GCCC 


GGAA 


GGAC 


GGAG 


GGAA 


GGCA 


GGAA 
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The oligonucleotide tags of the invention and their complements are 
conveniently synthesized on an automated DNA synthesizer, e.g. an Applied 
Biosystems, Inc. (Foster City, California) model 392 or 394 DNA/RNA Synthesizer, 
using standard chemistries, such as phosphoramidite chemistry, e.g. disclosed in the 
5 following references: Beaucage and Iyer, Tetrahedron, 48: 2223-231 1 (1992); Molko 
et al, U.S. patent 4,980,460; Koster et al, U.S. patent 4,725,677; Caruthers et al, U.S. 
patents 4,415,732; 4,458,066; and 4,973,679; and the like. Alternative chemistries, 
e.g. resulting in non-natural backbone groups, such as phosphorothioate, 
phosphoramidate, and the like, may also be employed provided that the resulting 

1 0 oligonucleotides are capable of specific hybridization. In some embodiments, tags 
may comprise naturally occurring nucleotides that permit processing or manipulation 
by enzymes, while the corresponding tag complements may comprise non-natural 
nucleotide analogs, such as peptide nucleic acids, or like compounds, that promote the 
formation of more stable duplexes during sorting. 

1 5 When microparticles are used as supports, repertoires of oligonucleotide tags 

and tag complements may be generated by subunit-wise synthesis via "split and mix" 
techniques, e.g. as disclosed in Shortle et al, International patent application 
PCT/US93/03418 or Lyttle et al, Biotechniques, 19: 274-280 (1995). Briefly, the 
basic unit of the synthesis is a subunit of the oligonucleotide tag. Preferably, 

20 phosphoramidite chemistry is used and 3* phosphoramidite oligonucleotides are 
prepared for each subunit in a minimally cross-hybridizing set, e.g. for the set first 
listed above, there would be eight 4-mer S-phosphoramidites. Synthesis proceeds as 
disclosed by Shortle et al or in direct analogy with the techniques employed to 
generate diverse oligonucleotide libraries using nucleosidic monomers, e.g. as 

25 disclosed in Telenius et al, Genomics, 13: 718-725 (1992); Welsh et al, Nucleic Acids 
Research, 19: 5275-5279 (1991); Grothues et al, Nucleic Acids Research, 21 : 1321- 
1322 (1993); Hartley, European patent application 90304496.4; Lam et al, Nature, 
354: 82-84 (1991); Zuckerman et al, Int. J. Pept. Protein Research, 40: 498-507 
(1992); and the like. Generall y, the se techniques sim pl y call for the a pplication of 

30 mixtures of the activated monomers to the growing oligonucleotide during the 

coupling steps. Preferably, oligonucleotide tags and tag complements are synthesized 
on a DNA synthesizer having a number of synthesis chambers which is greater than or 
equal to the number of different kinds of words used in the construction of the tags. 
That is, preferably there is a synthesis chamber corresponding to each type of word. 

35 In this embodiment, words are added nucleotide-by-nucleotide, such that if a word 
consists of five nucleotides there are five monomer couplings in each synthesis 
chamber. After a word is completely synthesized, the synthesis supports are removed 
from the chambers, mixed, and redistributed back to the chambers for the next cycle 
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of word addition. This latter embodiment takes advantage of the high coupling yields 
of monomer addition, e.g. in phosphoramidite chemistries. 

Double stranded forms of tags may be made by separately synthesizing the 
complementary strands followed by mixing under conditions that permit duplex 
5 formation. Alternatively, double stranded tags may be formed by first synthesizing a 
single stranded repertoire linkied to a known oligonucleotide sequence that serves as a 
primer binding site. The second strand is then synthesized by combining the single 
stranded repertoire with a primer and extending with a polymerase. This latter 
approach is described in Oliphant et al, Gene, 44: 177-183 (1986). Such duplex tags 

1 0 may then be inserted into cloning vectors along with target polynucleotides for sorting 
and manipulation of the target polynucleotide in accordance with the invention. 

When tag complements are employed that are made up of nucleotides that 
have enhanced binding characteristics, such as PNAs or oligonucleotide N3'— >P5' 
phosphoramidates, sorting can be implemented through the formation of D-loops 

1 5 between tags comprising natural nucleotides and their PNA or phosphoramidate 
complements, as an alternative to the "stripping" reaction employing the 3'— »5' 
exonuclease activity of a DNA polymerase to render a tag single stranded. 

Oligonucleotide tags of the invention may range in length from 12 to 60 
nucleotides or basepairs. Preferably, oligonucleotide tags range in length from 1 8 to 

20 40 nucleotides or basepairs. More preferably, oligonucleotide tags range in length 
from 25 to 40 nucleotides or basepairs. In terms of preferred and more preferred 
numbers of subunits, these ranges may be expressed as follows: 

Table IV 

25 Numbers of Subunits in Tags in Preferred Embodiments 



Monomers 

in Subunit Nucleotides in Oligonucleotide Tag 

(12-60) (18-40) (25-40) 



3 


4-20 subunits 


6-13 subunits 


8-13 subunits 


4 


3-15 subunits 


4-10 subunits 


6-10 subunits 


5 


2-12 subunits 


3-8 subunits 


5-8 subunits 


6 


2-10 subunits 


3-6 subunits 


4-6 subunits 



30 Most preferably, oligonucleotide tags are single stranded and specific hybridization 
occurs via Watson-Crick pairing with a tag complement. 
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Preferably, repertoires of single stranded oligonucleotide tags of the invention 
contain at least 1 00 members; more preferably, repertoires of such tags contain at 
least 1000 members; and most preferably, repertoires of such tags contain at least 
10,000 members. 

5 

Triplex Tags 

In embodiments where specific hybridization occurs via triplex formation, 
coding of tag sequences follows the same principles as for duplex-forming tags; 
however, there are further constraints on the selection of subunit sequences. 

1 0 Generally, third strand association via Hoogsteen type of binding is most stable along 
homopyrimidine-homopurine tracks in a double stranded target. Usually, base triplets 
form in T-A*T or C-G*C motifs (where indicates Watson-Crick pairing and "*" 
indicates Hoogsteen type of binding); however, other motifs are also possible. For 
example, Hoogsteen base pairing permits parallel and antiparallel orientations 

1 5 between the third strand (the Hoogsteen strand) and the purine-rich strand of the 

duplex to which the third strand binds, depending on conditions and the composition 
of the strands. There is extensive guidance in the literature for selecting appropriate 
sequences, orientation, conditions, nucleoside type (e.g. whether ribose or 
deoxyribose nucleosides are employed), base modifications (e.g. methylated cytosine, 

20 and the like) in order to maximize, or otherwise regulate, triplex stability as desired in 
particular embodiments, e.g. Roberts et al, Proc. Natl. Acad. Sci., 88: 9397-9401 
(1991); Roberts et al, Science, 258: 1463-1466 (1992); Roberts et al, Proc. Natl. 
Acad. Sci., 93: 4320-4325 (1996); Distefano et al, Proc. NatTAcadrSci., 90: 1 179- 
1 183 (1993); Mergny et al, Biochemistry, 30: 9791-9798 (1991); Cheng et al, J. Am. 

25 Chem. Soc, 1 14: 4465-4474 (1992); Beal and Dervan, Nucleic Acids Research, 20: 
2773-2776 (1992); Beal and Dervan, J. Am. Chem. Soc, 1 14: 4976-4982 (1992); 
Giovannangeli et al, Proc. Natl. Acad. Sci., 89: 8631-8635 (1992); Moser and Dervan, 
Science, 238: 645-650 (1987); McShan et al, J. Biol. Chem., 267:5712-5721 (1992); 

Yoon-etal^roc.-Nad..Ac^ al, Nucleic Acids 

30 Research, 20: 1777-1784 (1992); Thuong and Helene, Angew. Chem. Int. Ed. Engl. 
32: 666-690 (1993); Escude et al, Proc. Natl. Acad. Sci., 93: 4365-4369 (1996); and 
the like. Conditions for annealing single-stranded or duplex tags to their single- 
stranded or duplex complements are well known, e.g. Ji et al, Anal. Chem. 65: 1323- 
1328 (1993); Cantor et al, U.S. patent 5,482,836; and the like. Use of triplex tags has 

35 the advantage of not requiring a "stripping" reaction with polymerase to expose the 
tag for annealing to its complement. 

Preferably, oligonucleotide tags of the invention employing triplex 
hybridization are double stranded DNA and the corresponding tag complements are 



-15- 



WO 96/41011 



PCT/US96/09513 



single stranded. More preferably, 5-methylcytosine is used in place of cytosine in the 
tag complements in order to broaden the range of pH stability of the triplex formed 
between a tag and its complement. Preferred conditions for forming triplexes are 
fully disclosed in the above references. Briefly, hybridization takes place in 
5 concentrated salt solution, e.g. 1 .0 M NaCl, 1 .0 M potassium acetate, or the like, at 
pH below 5.5 ( or 6.5 if 5-methylcytosine is employed). Hybridization temperature 
depends on the length and composition of the tag; however, for an 1 8-20-mer tag of 
longer, hybridization at room temperature is adequate. Washes may be conducted 
with less concentrated salt solutions, e.g. 10 mM sodium acetate, 100 mM MgCl 2 , pH 

10 5.8, at room temperature. Tags may be eluted from their tag complements by 
incubation in a similar salt solution at pH 9.0. 

Minimally cross-hybridizing sets of oligonucleotide tags that form triplexes 
may be generated by the computer program of Appendix Ic, or similar programs. An 
exemplary set of double stranded 8-mer words are listed below in capital letters with 

1 5 the corresponding complements in small letters. Each such word differs from each of 
the other words in the set by three base pairs. 

Table V 

Exemplary Minimally Cross-Hvbridizing 
20 Set of DoubleStranded 8-mer Tags 



5' 


-AAGGAGAG 


5' 


-AAAGGGGA 


5' 


-AGAGAAGA 


5' 


-AGGGGGGG 


3' 


-TTCCTCTC 


3' 


-TTTCCCCT 


3' 


-TCTCTTCT 


3' 


-TCCCCCCC 


3' 


— ttcctctC 


3' 


-tttcccct 


3' 


-tctcttct 


3' 


-tccccccc 


5' 


-AAAAAAAA 


5' 


-AAGAGAGA 


5' 


-AGGAAAAG 


5' 


-GAAAGGAG 


3' 




3' 


-TTCTCTCT 


3' 


-TCCTTTTC 


3' 


-CTTTCCTC 


3' 


-tttttttt 


3' 


-ttctctct 


3' 


-tccttttc 


3' 


-ctttcctc 


5' 


-AAAAAGGG 


5' 


-AGAAGAGG 


5 r 


-AGGAAGGA 


5' 


-GAAGAAGG 


3' 


-TTTTTCCC 


3' 


-TCTTCTCC 


3' 


-TCCTTCCT 


3' 


-CTTCTTCC 


3' 


-tttttccc 


3' 


-tcttctcc 


3' 


-tccttcct 


3' 


-cttcttcc 


5' 


-AAAGGAAG 


5' 


-AGAAGGAA 


5' 


-AGGGGAAA 


5' 


-GAAGAGAA 
-CTTCTCTT" 


3' 


-T-T-TCCTTC 

-tttccttc 


3' 


-tcttcctt 


3' 


-tccccttt 


3' 


-cttctctt ' 



Table VI 

25 



Repertoire Size of Various Double Stranded Tags 
That Form Triplexes with Their Tag Complements 
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Oltgonucleotid 
c 

Word 
Length 


Nucleotide ■ 
Difference 
between 
Oligonucleotides 
of Minimally 
Cross- 
Hybridizing Set 


Maximal Size 
of Minimally 

Cross- 
Hybridizing 
Set 


Size of 
Repertoire 
with Four 

Words 


Size of 
Repertoire with 
Five Words 


4 


2 


8 


4096 


3.2 x 10 q 


6 


3 


8 


4096 


3.2x10" 


8 


3 


16 


6.5 x 10* 


1.05 x 10° 


10 


5 


8 


4096 




15 


5 


92 






20 


6 


765 






20 


8 


92 






20 


10 


22 







Preferably, repertoires of double stranded oligonucleotide tags of the invention 
contain at least 10 members; more preferably, repertoires of such tags contain at least 
100 members. Preferably, words are between 4 and 8 nucleotides in length for 
5 combinatorial^ synthesized double stranded oligonucletide tags, and oligonucleotide 
tags are between 12 and 60 base pairs in length. More preferably, such tags are 
between 1 8 and 40 base pairs in length. 

Solid Phase Supports 

10 Solid phase supports for use with the invention may have a wide variety of 

forms, including microparticles, beads, and membranes, slides, plates, micromachined 
chips, and the like. Likewise, solid phasejsupports of the invention may comprise a 
wide variety of compositions, including glass, plastic, silicon, alkanethiolate- 
derivatized gold, cellulose, low cross-linked and high cross-linked polystyrene, silica 

1 5 gel, polyamide, and the like. Preferably, either a population of discrete particles are 
employed such that each has a uniform coating, or population, of complementary 
sequences of the same tag (and no other), or a single or a few supports are employed 
with spatially discrete regions each containing a uniform coating, or population, of 
complementary sequences to the same tag (and no other). In the latter embodiment, 

20 the area of the regions may vary according to particular applications; usually, the 
regions range in area from several nm 2 , e.g. 3-5, to several hundred nm 2 , e.g. 100- 
500. Preferably, such regions are spatially discrete so that signals generated by 
events, e.g. fluorescent emissions, at adjacent regions can be resolved by the detection 
system being employed. In some applications, it may be desirable to have regions 

25 with uniform coatings of more than one tag complement, e.g. for simultaneous 

sequence analysis, or for bringing separately tagged molecules into close proximity. 
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Tag complements may be used with the solid phase support that they are 
synthesized on, or they may be separately synthesized and attached to a solid phase 
support for use, e.g. as disclosed by Lund et al, Nucleic Acids Research, 16: 10861- 
10880 (1988); Albretsen et al, Anal. Biochem., 189: 40-50 (1990); Wolf et al, Nucleic 
5 Acids Research, 1 5: 291 1-2926 (1987); or Ghosh et al/Nucleic Acids Research, 15: 
5353-5372 (1987). Preferably, tag complements are synthesized on and used with the 
same solid phase support, which may comprise a variety of forms and include a 
variety of linking moieties. Such supports may comprise microparticles or arrays, or 
matrices, of regions where uniform populations of tag complements are synthesized. 

10 A wide variety of microparticle supports may be used with the invention, including 
microparticles made of controlled pore glass (CPG), highly cross-linked polystyrene, 
acrylic copolymers, cellulose, nylon, dextran, latex, polyacrolein, and the like, 
disclosed in the following exemplary references: Meth. Enzymol., Section A, pages 
1 1-147, vol. 44 (Academic Press, New York, 1976); U.S. patents 4,678,814; 

1 5 4,413,070; and 4,046;720; and Pon, Chapter 19, in Agrawal, editor, Methods in 
Molecular Biology, Vol. 20, (Humana Press, Totowa, NJ, 1993). Microparticle 
supports further include commercially available nucleoside-derivatized CPG and 
polystyrene beads (e.g. available from Applied Biosystems, Foster City, CA); 
derivatized magnetic beads; polystyrene grafted with polyethylene glycol (e.g., 

20 TentaGelTM, Rapp Polymere, Tubingen Germany); and the like. Selection of the 
support characteristics, such as material, porosity, size, shape, and the like, and the 
type of linking moiety employed depends on the conditions under which the tags are 
used. For example, in applications involving successive processing with enzymes, 
supports and linkers that minimize steric hindrance of the enzymes and that facilitate 

25 access to substrate are preferred. Other important factors to be considered in selecting 
the most appropriate microparticle support include size uniformity, efficiency as a 
synthesis support, degree to which surface area known, and optical properties, e.g. as 
explain more fully below, clear smooth beads provide instrumentational advantages 

when handling large numbers of beads on a surfa ce. 

30 Exemplary linking moieties for attaching and/or synthesizing tags on 

microparticle surfaces are disclosed in Pon et al, Biotechniques, 6:768-775 (1988); 
Webb, U.S. patent 4,659,774; Barany et al, International patent application 
PCT/US91/06103; Brown et al, J. Chem. Soc. Commun., 1989: 891-893; Damha et 
al, Nucleic Acids Research, 18: 3813-3821 (1990); Beattie et al, Clinical Chemistry, 

35 39: 719-722 (1993); Maskos and Southern, Nucleic Acids Research, 20: 1679-1684 
(1992); and the like. 

As mentioned above, tag complements may also be synthesized on a single 
(or a few) solid phase support to form an array of regions uniformly coated with tag 
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complements. That is, within each region in such an array the same tag complement 
is synthesized. Techniques for synthesizing such arrays are disclosed in McGall et al, 
International application PCT/US93/03767; Pease et al, Proc. Natl. Acad. Sci.,-91: 
5022-5026 (1994); Southern and Maskos, International application 
5 PCT/GB89/01 1 14; Maskos and Southern (cited above); -Southern et al, Genomics, 13: 
1008-1017 (1992); and Maskos and Southern, Nucleic Acids Research, 21 : 4663- 
4669(1993). 

Preferably, the invention is implemented with microparticles or beads 
uniformly coated with complements of the same tag sequence. Microparticle supports 

1 0 and methods of covalently or noncovalently linking oligonucleotides to their surfaces 
are well known, as exemplified by the following references: Beaucage and Iyer (cited 
above); Gait, editor, Oligonucleotide Synthesis: A Practical Approach (IRL Press, 
Oxford, 1984); and the references cited above. Generally, the size and shape of a 
microparticle is not critical; however, microparticles in the size range of a few, e.g. 1- 

15 2, to several hundred, e.g. 200-1000 urn diameter are preferable, as they facilitate the 
construction and manipulation of large repertoires of oligonucleotide tags with 
minimal reagent and sample usage. 

In some preferred applications, commercially available controlled-pore glass 
(CPG) or polystyrene supports are employed as solid phase supports in the invention. 

20 Such supports come available with base-labile linkers and initial nucleosides attached, 
e.g. Applied Biosystems (Foster City, CA). Preferably, microparticles having pore 
size between 500 and 1000 angstroms are employed. 

In other preferred applications, non-porous microparticles are employed for 
their optical properties, which may be advantageously used when tracking large 

25 numbers of microparticles on planar supports, such as a microscope slide. 

Particularly preferred non-porous microparticles are the glycidal methacrylate (GMA) 
beads available from Bangs Laboratories (Carmel, IN). Such microparticles are 
useful in a variety of sizes and derivatized with a variety of linkage groups for 

synthesizing tag* nr ta g complements. Preferably^for massively-parallel 

30 manipulations of tagged microparticles, 5 urn diameter GMA beads are employed. 

Attaching Tags to Polynucleotides 
For Sorting onto Solid Phase Supports 
An important aspect of the invention is the sorting and attachment of a 
35 populations of polynucleotides, e.g. from a cDNA library, to microparticles or to 
separate regions on a solid phase support such that each microparticle or region has 
substantially only one kind of polynucleotide attached. This objective is 
accomplished by insuring that substantially all different polynucleotides have 
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different tags attached. This condition, in turn, is brought about by taking a sample of 
the full ensemble of tag-polynucleotide conjugates for analysis. (It is acceptable that 
identical polynucleotides have different tags, as it merely results in the same 
polynucleotide being operated on or analyzed twice in two different locations.) Such 
5 sampling can be carried out either overtly—for example, by taking a small volume 
from a larger mixture-after the tags have been attached to the polynucleotides, it can 
be carried out inherently as a secondary effect of the techniques used to process the 
polynucleotides and tags, or sampling can be carried out both overtly and as an 
inherent part of processing steps. 

10 Preferably, in constructing a cDNA library where substantially all different 

cDNAs have different tags, a tag repertoire is employed whose complexity, or number 
of distinct tags, greatly exceeds the total number of mRNAs extracted from a cell or 
tissue sample. Preferably, the complexity of the tag repertoire is at least 10 times that 
of the polynucleotide population; and more preferably, the complexity of the tag 

1 5 repertoire is at least 1 00 times that of the polynucleotide population. Below, a 
protocol is disclosed for cDNA library construction using a primer mixture that 
contains a full repertoire of exemplary 9-word tags. Such a mixture of tag-containing 
primers has a complexity of 8 9 , or about 1 .34 x 1 0 8 . As indicated by Winslow et al, 
Nucleic Acids Research, 19: 3251-3253 (1991), mRNA for library construction can 

20 be extracted from as few as 10-100 mammalian cells. Since a single mammalian cell 
contains about 5 x 10 5 copies of mRNA molecules of about 3.4 x 10 4 different kinds, 
by standard techniques one can isolate the mRNA from about 100 cells, or 
(theoretically) about 5 x 10 7 mRNA molecules. Comparing this number to the 
complexity of the primer mixture shows that without any additional steps, and even 

25 assuming that mRNAs are converted into cDNAs with perfect efficiency (1 % 

efficiency or less is more accurate), the cDNA library construction protocol results in 
a population containing no more than 37% of the total number of different tags. That 
is, without any overt sampling step at all, the protocol inherently generates a sample 

that comprises-37-%, or less, of the tag repertoire. The probability of obtaining a 

30 double under these conditions is about 5%, which is within the preferred range. With 
mRNA from 10 cells, the fraction of the tag repertoire sampled is reduced to only 
3.7%, even assuming that all the processing steps take place at 100% efficiency. In 
fact, the efficiencies of the processing steps for constructing cDNA libraries are very 
low, a "rule of thumb" being that good library should contain about 10 8 cDNA clones 

35 from mRNA extracted from 1 0 6 mammalian cells. 

Use of larger amounts of mRNA in the above protocol, or for larger amounts 
of polynucleotides in general, where the number of such molecules exceeds the 
complexity of the tag repertoire, a tag-polynucleotide conjugate mixture potentially 
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contains every possible pairing of tags and types of mRNA or polynucleotide. In such 
cases, overt sampling may be implemented by removing a sample volume after a 
serial dilution of the starting mixture of tag-polynucleotide conjugates. The amount 
of dilution required depends on the amount of starting material and the efficiencies of 
5 the processing steps, which are readily estimated. 

If mRNA were extracted from 10 6 cells (which would correspond to about 0.5 
ug of poly(A) + RNA), and if primers were present in about 10-100 fold concentration 
excess-as is called for in a typical protocol, e.g. Sambrook et al, Molecular Cloning, 
Second Edition, page 8.61 [10 |iL 1.8 kb mRNA at 1 mg/mL equals about 1.68 x 10"" 
1 0 moles and 1 0 uL 1 8-mer primer at 1 mg/mL equals about 1 .68 x 1 0* 9 moles], then the 
total number of tag-polynucleotide conjugates in a cDN A library would simply be 
equal to or less than the starting number of mRNAs, or about 5 x 10 n vectors 
containing tag-polynucleotide conjugates-again this assumes that each step in cDNA 
construction-first strand synthesis, second strand synthesis, ligation into a vector- 
1 5 occurs with perfect efficiency, which is a very conservative estimate. The actual 
number is significantly less. 

If a sample of n tag-polynucleotide conjugates are randomly drawn from a 
reaction mixture-as could be effected by taking a sample volume, the probability of 
drawing conjugates having the same tag is described by the Poisson distribution, 
20 P(r)=e"*XX)7r, where r is the number of conjugates having the same tag and X=np, 
where p is the probability of a given tag being selected. If n=10 6 and p=l/(1.34 x 
10 8 ), then X=.00746 and P(2)=2.76 x 10" 5 . Thus, a sample of one million molecules 
gives rise to an expected number of doubles well within the preferred range. Such a 
sample is readily obtained as follows: Assume that the 5 x 10 n mRNAs are perfectly 
25 converted into 5 x 10 u vectors with tag-cDNA conjugates as inserts and that the 5 x 
1 0 1 1 vectors are in a reaction solution having a volume of 1 00 ul. Four 1 0-fold serial 
dilutions may be carried out by transferring 10 ul from the original solution into a 
vessel containing 90 ul of an appropriate buffer, such as TE. This process may be 

repeated forthree-atiditior^ — 

30 vector molecules per ul. A 2 ul aliquot from this solution yields 1 0 6 vectors 

containing tag-cDNA conjugates as inserts. This sample is then amplified by straight 
forward transformation of a competent host cell followed by culturing. 

Of course, as mentioned above, no step in the above process proceeds with 
perfect efficiency. In particular, when vectors are employed to amplify a sample of 
35 tag-polynucleotide conjugates, the step of transforming a host is very inefficient. 
Usually, no more than 1% of the vectors are taken up by the host and replicated. 
Thus, for such a method of amplification, even fewer dilutions would be required to 
obtain a sample of 10 6 conjugates. 
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A repertoire of oligonucleotide tags can be conjugated to a population of 
polynucleotides in a number of ways, including direct enzymatic ligation, 
amplification, e.g. via PCR, using primers containing the tag sequences, and the like. 
The initial ligating step produces a very large population of tag-polynucleotide 
5 conjugates such that a single tag is generally attached to-many different 

polynucleotides. However, as noted above, by taking a sufficiently small sample of 
the conjugates, the probability of obtaining "doubles," i.e. the same tag on two 
different polynucleotides, can be made negligible. Generally, the larger the sample 
the greater the probability of obtaining a double. Tlius, a design trade-off exists 

1 0 between selecting a large sample of tag-polynucleotide conjugates-which, for 
example, ensures adequate coverage of a target polynucleotide in a shotgun 
sequencing operation or adequate representation of a rapidly changing mRNA pool, 
and selecting a small sample which ensures that a minimal number of doubles will be 
present. In most embodiments, the presence of doubles merely adds an additional 

1 5 source of noise or, in the case of sequencing, a minor complication in scanning and 
signal processing, as microparticles giving multiple fluorescent signals can simply be 
ignored. 

As used herein, the term "substantially all" in reference to attaching tags to 
molecules, especially polynucleotides, is meant to reflect the statistical nature of the 

20 sampling procedure employed to obtain a population of tag-molecule conjugates 
essentially free of doubles. The meaning of substantially all in terms of actual 
percentages of tag-molecule conjugates depends on how the tags are being employed. 
Preferably, for nucleic acid sequencing, substantially all means that at least eighty 
percent of the polynucleotides have unique tags attached. More preferably, it means 

25 that at least ninety percent of the polynucleotides have unique tags attached. Still 

more preferably, it means that at least ninety-five percent of the polynucleotides have 
unique tags attached. And, most preferably, it means that at least ninety-nine percent 
of the polynucleotides have unique tags attached. 

Preferably , when the populatio n of polynucleotides consists of messenger 

30 RNA (mRNA), oligonucleotides tags may be attached by reverse transcribing the 
mRNA with a set of primers preferably containing complements of tag sequences. 
An exemplary set of such primers could have the following sequence: 

5 1 -mRNA- [A] n -3 1 
35 [T ] ! 9 GG [ W, W, W, C] qAC CAGCTGA TC-5 1 -biotin 
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where "[W,W,W,C]9 M represents the sequence of an oligonucleotide tag of nine 
subunits of four nucleotides each and n [W,W,W,C]" represents the subunit sequences 
listed above, i.e. "W" represents T or A. The underlined sequences identify an 
optional restriction endonuclease site that can be used to release the polynucleotide 
5 from attachment to a solid phase support via the biotin, if one is employed. For the 
above primer, the complement attached to a microparticle could have the form: 

5 ' - [G, W, W, W] 9 TGG-linker-microparticle 

1 0 After reverse transcription, the mRNA is removed, e.g. by RNase H digestion, 

and the second strand of the cDNA is synthesized using, for example, a primer of the 
following form: 

5 ' -NRRGATCYNNN-3 1 

15 

where N is any one of A, T, G, or C; R is a purine-containing nucleotide, and Y is a 
pyrimidineKxmtaining nucleotide. This particular primer creates a Bst Yl restriction 
site in the resulting double stranded DN A which, together with the Sal I site, 
facilitates cloning into a vector with, for example, Bam HI and Xho I sites. After Bst 
20 Yl and Sal I digestion, the exemplary conjugate would have the form: 

S'-RCGACCAlCW.W.WlgGGCTJxg- cDNA -NNNR 

GGT[G,W,W,W] 9 CC[A] 19 - rDNA -NNNYCTAG-5' 

25 The polynucleotide-tag conjugates may then be manipulated using standard molecular 
biology techniques. For example, the above conjugate-which is actually a mixture- 
may be inserted into commercially available cloning vectors, e.g. Stratagene Cloning 
System (La Jolla, CA); transfected into a host, such as a commercially available host 
bacteria;~which-is-thenxulture^ 

30 vectors may then be isolated using standard techniques, e.g. Sambrook et al, 

Molecular Cloning, Second Edition (Cold Spring Harbor Laboratory, New York, 
1989). Alternatively, appropriate adaptors and primers may be employed so that the 
conjugate population can be increased by PCR. 

Preferably, when the ligase-based method of sequencing is employed, the Bst 

35 Yl and Sal I digested fragments are cloned into a Bam HWXho I-digested vector 
having the following single-copy restriction sites: 
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5 1 -GA GGATG CCTTTAT GGATCCACTCGAGA TCCCAATCCA-3 1 
Fokl BamHI Xhol 

This adds the Fok I site which will allow initiation of the sequencing process 
5 discussed more fully below. 

Tags can be conjugated to cDNAs of existing libraries by standard cloning 
methods. cDNAs are excised from their existing vector, isolated, and then ligated into 
a vector containing a repertoire of tags. Preferably, the tag-containing vector is 
linearized by cleaving with two restriction enzymes so that the excised cDNAs can be 

10 ligated in a predetermined orientation. The concentration of the linearized tag- 
containing vector is in substantial excess over that of the cDNA inserts so that 
ligation provides an inherent sampling of tags. 

A general method for exposing the single stranded tag after amplification 
involves digesting a target polynucleotide-containing conjugate with the 5'— >3* 

1 5 exonuclease activity of T4 DNA polymerase, or a like enzyme. When used in the 
presence of a single deoxynucleoside triphosphate, such a polymerase will cleave 
nucleotides from 3* recessed ends present on the non-template strand of a double 
stranded fragment until a complement of the single deoxynucleoside triphosphate is 
reached on the template strand. When such a nucleotide is reached the 5 — >3* 

20 digestion effectively ceases, as the polymerase's extension activity adds nucleotides at 
a higher rate than the excision activity removes nucleotides. Consequently, single 
stranded tags constructed with three nucleotides are readily prepared for loading onto 
solid phase supports. 

The technique may also be used to preferentially methylate interior Fok I sites 

25 of a target polynucleotide while leaving a single Fok I site at the terminus of the 

polynucleotide unmethylated. First, the terminal Fok I site is rendered single stranded 
using a polymerase with deoxycytidine triphosphate. The double stranded portion of 
the fragment is then methylated, after which the single stranded terminus is filled in 

with a DNA polymerase in the presence of all four nucleoside triphosphates, ther eby 

30 regenerating the Fok I site. Clearly, this procedure can be generalized to 
endonucleases other than Fok I. 

After the oligonucleotide tags are prepared for specific hybridization, e.g. by 
rendering them single stranded as described above, the polynucleotides are mixed 
with microparticles containing the complementary sequences of the tags under 

35 conditions that favor the formation of perfectly matched duplexes between the tags 

and their complements. There is extensive guidance in the literature for creating these 
conditions. Exemplary references providing such guidance include Wetmur, Critical 
Reviews in Biochemistry and Molecular Biology, 26: 227-259 (1991); Sambrook et 
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al. Molecular Cloning: A Laboratory Manual, 2nd Editic. Cold Spring Harbor 
Laboratory, New York, 1989); and the like. Preferably, uv. hybridization conditions 
are sufficiently stringent so that only perfectly matched sequences form stable 
duplexes. Under such conditions the polynucleotides specifically hybridized through 
5 their tags may be ligated to the complementary sequences attached to the 

microparticles. Finally, the microparticles are washed to remove polynucleotides with 
unligated and/or mismatched tags. 

When CPG microparticles conventionally employed as synthesis supports are 
used, the density of tag complements on the microparticle surface is typically greater 
1 0 than that necessary for some sequencing operations. That is, in sequencing 

approaches that require successive treatment of the attached polynucleotides with a 
variety of enzymes, densely spaced polynucleotides may tend to inhibit access of the 
relatively bulky enzymes to the polynucleotides. In such cases, the polynucleotides 
are preferably mixed with the microparticles so that tag complements are present in 
1 5 significant excess, e.g. from 10:1 to 1 00: 1 , or greater, over the polynucleotides. This 
ensures that the density of polynucleotides on the microparticle surface will not be so 
high as to inhibit enzyme access. Preferably, the average inter-polynucleotide spacing 
on the microparticle surface is on the order of 30-100 nm. Guidance in selecting 
ratios for standard CPG supports and Ballotini beads (a type of solid glass support) is . 
20 found in Maskos and Southern, Nucleic Acids Research, 20: 1 679- 1 684 (1 992). 

Preferably, for sequencing applications, standard CPG beads of diameter in the range 
of 20-50 ^im are loaded with about 10 5 polynucleotides, and GMA beads of diameter 
in the range of 5-10 nm are loaded with a few tens of thousand polynucleotide, e.g. 4 
xK^toexlO 4 . 

25 In the preferred embodiment, tag complements are synthesized on 

microparticles combinatorially; thus, at the end of the synthesis, one obtains a 
complex mixture of microparticles from which a sample is taken for loading tagged 
polynucleotides. The size of the sample of microparticles will depend on several 



faetorsrineludmg-the-size-o£the 

30 apparatus for used for observing loaded microparticles-e.g. its capacity, the tolerance 
for multiple copies of microparticles with the same tag complement (i.e. "bead 
doubles"), and the like. The following table provide guidance regarding 
microparticle sample size, microparticle diameter, and the approximate physical 
dimensions of a packed array of microparticles of various diameters. 




35 
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Micropaiticle diameter 



5 fim 



10 jim 20 40 um 



Max. no. 

polynucleotides loaded 
at 1 per 10 5 sq. 
angstrom 

Approx. area of 
monolayer of 10 6 
microparticles 



3xl0 5 1.26 xlO 6 5xl0 6 



.45 x .45 cm 1 x 1 cm 2 x 2 cm 4 x 4 cm 



The probability that the sample of microparticles contains a given tag complement or 
is present in multiple copies is described by the Poisson distribution, as indicated in 
the following table. 



Table VII 



Number of 
microparticles in 
sample (as fraction 
of repertoire size), 



m 



1.000 
.693 
.405 
.285 
.223 
.105 
.010 



Fraction of 
repertoire of tag 
complements 
present in 
sample, 



0.63 
0.50 
0.33 
0.25 
0.20 
0.10 
0.01 



Fraction of 
microparticles in 
sample with unique 
tag complement 
attached, 



m(e' m )/2 



0.37 
0.35 
0.27 
0.21 
0.18 
0.09 
0.01 



Fraction of 
microparticles in 
sample carrying 
same tag 
complement as one 
other micropaiticle 

in sample 
("bead doubles"), 
mHe m )/2 



0.18 
0.12 
0.05 
0.03 
0.02 
0.005 



High Specificity Sort ing and Panning 
The kinetics of sorting depends on the rate of hybridization of oligonucleotide 

10 tags to their tag complements which, in turn, depends on the complexity of the tags in 
the hybridization reaction. Thus, a trade off exists between sorting rate and tag 
complexity, such that an increase in sorting rate may be achieved at the cost of 
reducing the complexity of the tags involved in the hybridization reaction. As 
explained below, the effects of this trade off may be ameliorated by "panning." 

1 5 Specificity of the hybridizations may be increased by taking a sufficiently 

small sample so that both a high percentage of tags in the sample are unique and the 
nearest neighbors of substantially all the tags in a sample differ by at least two words. 
This latter condition may be met by taking a sample that contains a number of tag- 
polynucleotide conjugates that is about 0.1 percent or less of the size of the repertoire 
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being employed. For example, if tags are constructed with eight words selected from 
Table II, a repertoire of 8 8 , or about 1 .67 x 1 0 7 , tags and tag complements are 
produced. In a library of tag-cDNA conjugates as described above, a 0.1 percent 
sample means that about 16,700 different tags are present. If this were loaded directly 
5 onto a repertoire-equivalent of microparticles, or in this example a sample of 1 .67 x 
10 7 microparticles, then only a sparse subset of the sampled microparticles would be 
loaded. The density of loaded microparticles can be increase-for example, for more 
efficient sequencing~by undertaking a "panning" step in which the sampled tag- 
cDNA conjugates are used to separate loaded microparticles from unloaded 

1 0 microparticles. Thus, in the example above, even though a "0.1 percent" sample 
contains only 16,700 cDNAs, the sampling and panning steps may be repeated until 
as many loaded microparticles as desired are accumulated. 

A panning step may be implemented by providing a sample of tag-cDNA 
conjugates each of which contains a capture moiety at an end opposite, or distal to, 

1 5 the oligonucleotide tag. Preferably, the capture moiety, is of a type which can be 
released from the tag-cDNA conjugates, so that the tag-cDNA conjugates can be 
sequenced with a single-base sequencing method. Such moieties may comprise 
biotin, digoxigenin, or like ligands, a triplex binding region, or the like. Preferably, 
such a capture moiety comprises a biotin component. Biotin may be attached to tag- 

20 cDN A conjugates by a number of standard techniques. If appropriate adapters 

containing PCR primer binding sites are attached to tag-cDNA conjugates, biotin may 
be attached by using a biotinylated primer in an amplification after sampling. 
Alternatively, if the tag-cDNA conjugates are inserts of cloning vectors, biotin may be 
attached after excising the tag-cDNA conjugates by digestion with an appropriate 

25 restriction enzyme followed by isolation and filling in a protruding strand distal to the 
tags with a DNA polymerase in the presence of biotinylated uridine triphosphate. 

After a tag-cDNA conjugate is captured, it may be released from the biotin 
moiety in a number of ways, such as by a chemical linkage that is cleaved by 

reduction > -e.g.-Herman-et al,-Anal^Bioch em ., 1 56: 4^55-(-L986) T -or-that-isxleaved 

30 photochemically, e.g. Olejnik et al, Nucleic Acids Research, 24: 361-366 (1996), or 
that is cleaved enzymatically by introducing a restriction site in the PCR primer. The 
latter embodiment can be exemplified by considering the library of tag-polynucleotide 
conjugates described above: 

35 5 ' -RCGACCA [C, W, W, W] 9GG [T] 19- cDNA -NNNR 

GGT[G,W,W,W] 9 CC[A) 19 - rDNA -NNNYCTAG-5' 
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The following adapters may be ligated to the ends of these fragments to permit 
amplification by PCR: 

5 • - xxxxxxxxxxxxxxxxxxxx 

5 XXXXXXXXXXXXXXXXXXXXYGAT 

Right Adapter 

GATCZZACTAGTZZZZZZZZZZZZ-3 ' 
10 ZZTGATCAZZZZZZZZZZZZ 

Left Adapter 
ZZTGATCAZZZZZZZZZZZZ-5 1 -biotin 

15 

Left Primer 

where "ACTAGT" is a Spe I recognition site (which leaves a staggered cleavage 
ready for single base sequencing), and the X's and Z's are nucleotides selected so that 
20 the annealing and dissociation temperatures of the respective primers are 

approximately the same. After ligation of the adapters and amplification by PCR 
using the biotinylated primer, the tags of the conjugates are rendered single stranded 
by the exonuclease activity of T4 DNA polymerase and conjugates are combined with 
a sample of microparticles, e.g. a repertoire equivalent, with tag complements 
25 attached. After annealing under stringent conditions (to minimize mis-attachment of 
tags), the conjugates are preferably ligated to their tag complements and the loaded 
microparticles are separated from the unloaded microparticles by capture with 
avidinated magnetic beads, or like capture technique. 

Returning to the example, this process results in the accumulation of about 
30 10,500 (=16,700 x .63) loaded microparticles with different tags, which may be 
released fronvme magnetic beads4>y^ 
40-50 times with new samples of microparticles and tag-cDNA conjugates, 4-5 x 10 5 
cDNAs can be accumulated by pooling the released microparticles. The pooled 
microparticles may then be simultaneously sequenced by a single-base sequencing 
35 technique. 

Determining how many times to repeat the sampling and panning steps-or 
more generally, deterrnining how many cDNAs to analyze, depends on one's 
objective. If the objective is to monitor the changes in abundance of relatively 
common sequences, e.g. making up 5% or more of a population, then relatively small 
40 samples, i.e. a small fraction of the total population size, may allow statistically 
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significant estimates of relative abundances. On the other hand, if one seeks to 
monitor the abundances of rare sequences, e.g. making up 0.1% or less of a 
population, then large samples are required. Generally, there is a direct relationship 
between sample size and the reliability of the estimates of relative abundances based 
5 on the sample. There is extensive guidance in the literature on determining 

appropriate sample sizes for making reliable statistical estimates, e.g. Koller et al, 
Nucleic Acids Research, 23:185-191 (1994); Good, Biometrika, 40: 16-264 (1953); 
Bunge et al, J. Am. Stat. Assoc., 88: 364-373 (1993); and the like. Preferably, for 
monitoring changes in gene expression based on the analysis of a series of cDNA 

10 libraries containing 10 5 to 10 8 independent clones of 3.0-3.5 x 10 4 different 

sequences, a sample of at least 10 4 sequences are accumulated for analysis of each 
library. More preferably, a sample of at least 10 5 sequences are accumulated for the 
analysis of each library; and most preferably, a sample of at least 5 x 10 5 sequences 
are accumulated for the analysis of each library. Alternatively, the number of 

1 5 sequences sampled is preferably sufficient to estimate the relative abundance of a 
sequence present at a frequency within the range of 0.1% to 5% with a 95% 
confidence limit no larger than 0.1% of the population size. 

Single Base DNA Sequencing 

20 The present invention can be employed with conventional methods of DNA 

sequencing, e.g. as disclosed by Hultman et al, Nucleic Acids Research, 17: 4937- 
4946 (1989). However, for parallel, or simultaneous, sequencing of multiple 
polynucleotides, a DNA sequencing niethodology is preferred that requires neither 
electrophoretic separation of closely sized DNA fragments nor analysis of cleaved 

25 nucleotides by a separate analytical procedure, as in peptide sequencing. Preferably, 
the methodology permits the stepwise identification of nucleotides, usually one at a 
time, in a sequence through successive cycles of treatment and detection. Such 
methodologies are referred to herein as "single base" sequencing methods. Single 
base-approaches-ar^disGlosed-in-the^ollowing-references^Cheesern^ p atent — 

30 5,302,509; Tsien et al, International application WO 91/06678; Rosenthal et al, 
International application WO 93/21340; Canard et al, Gene, 148: 1-6 (1994); and 
Metzker et al, Nucleic Acids Research, 22: 4259-4267 (1994). 

A "single base" method of DNA sequencing which is suitable for use with the 
present invention and which requires no electrophoretic separation of DNA fragments 

35 is described in International application PCT/US95/03678. Briefly, the method 
comprises the following steps: (a) ligating a probe to an end of the polynucleotide 
having a protruding strand to form a ligated complex, the probe having a 
complementary protruding strand to that of the polynucleotide and the probe having a 
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nuclease recognition site; (b) removing unligated probe from the ligated complex; (c) 
identifying one or more nucleotides in the protruding strand of the polynucleotide by 
the identity of the ligated probe; (d) cleaving the ligated complex with a nuclease; and 
(e) repeating steps (a) through (d) until the nucleotide sequence of the polynucleotide, 
5 or a portion thereof, is determined. 

A single signal generating moiety, such as a single fluorescent dye, may be 
employed when sequencing several different target polynucleotides attached to 
different spatially addressable solid phase supports, such as fixed microparticles, in a 
parallel sequencing operation. This may be accomplished by providing four sets of 
1 0 probes that are applied sequentially to the plurality of target polynucleotides on the 
different microparticles. An exemplary set of such probes are shown below: 



Set] 



Set 2 



Set 3 



Set 4 



ANNNN. . . NN 

N. . .NNTT. . .T* 

dCNNNN. . .NN 

N. . .NNTT. . .T 

dGNNNN. . .NN 

N. . .NNTT. . .T 

dTNNNN . . . NN 

N NNTT. . .T 



dANNNN. . .NN dANNNN. . . NN dANNNN . . . NN 

d N...NNTT...T N...NNTT...T N. . .NNTT. . .T 

CNNNN . . .NN dCNNNN. . .NN dCNNNN . . . NN 

N. . .NNTT. . .T* N. . .NNTT. . ,T N. . .NNTT. . .T 

dGNNNN... NN GNNNN . . . NN dGNNNN. . .NN 

N . . . NNTT T N . . . NNTT . . . T* N . . . NNTT . . . T 

dTNNNN... NN dTNNNN. . .NN TNNNN. . .NN 

N...NNTT...T N...NNTT...T N. . .NNTT. . . T* 



where each of the listed probes represents a mixture of 4 3 =64 oligonucleotides such 
that the identity of the 3' terminal nucleotide of the top strand is fixed and the other 
positions in the protruding strand are filled by every 3-mer permutation of nucleotides, 
or complexity reducing analogs. The listed probes are also shown with a single 
stranded poly-T tail with a signal generating moiety attached to the terminal thymidine, 
shown as "T* H . The "d" on the unlabeled probes designates a ligation-blocking moiety 
or absense ot 3'-hydroxyl, which prevents unlabeled probes from being ligated. 
Preferably, such 3-terminal nucleotides are dideoxynucleotides. In this embodiment, 
the probes of set lare first applied to the plurality of target polynucleotides and treated 
with a ligase so that target polynucleotides having a thymidine complementary to the 3' 
25 terminal adenosine of the labeled probes are ligated. The unlabeled probes are 

simultaneously applied to minimize inappropriate ligations. The locations of the target 
polynucleotides that form ligated complexes with probes terminating in "A" are 
identified by the signal generated by the label carried on the probe. After washing and 
cleavage, the probes of set 2 are applied. In this case, target polynucleotides forming 



15 



20 



-30- 



WO 96/41011 



PCT/US96/09513 



ligated complexes with probes terminating in "C" are identified by location. Similarly, 
the probes of sets 3 and 4 are applied and locations of positive signals identified. This 
process of sequentially applying the four sets of probes continues until the desired 
number of nucleotides are identified on the target polynucleotides. Clearly, one of 
5 ordinary skill could construct similar sets of probes that could have many variations, 
such as having protruding strands of different lengths, different moieties to block 
ligation of unlabeled probes, different means for labeling probes, and the like. 

A pparatus for Observing Enzymatic Processes and/or 

10 Binding Events at Microparticle Surfaces 

An objective of the invention is to sort identical molecules, particularly 
polynucleotides, onto the surfaces of microparticles by the specific hybridization of 
tags and their complements. Once such sorting has taken place, the presence of the 
molecules or operations performed on them can be detected in a number of ways 

1 5 depending on the nature of thelagged molecule, whether microparticles are detected 
separately or in "batches," whether repeated measurements are desired, and the like. 
Typically, the sorted molecules are exposed to ligands for binding, e.g. in drug 
development, or are subjected chemical or enzymatic processes, e.g. in polynucleotide 
sequencing. In both of these uses it is often desirable to simultaneously observe 

20 signals corresponding to such events or processes on large numbers of microparticles. 
Microparticles carrying sorted molecules (referred to herein as "loaded" 
microparticles) lend themselves to such large scale parallel operations, e.g. as 
demonstrated by Lam et al (cited above). ~ 

Preferably, whenever light-generating signals, e.g. chemiluminescent, 

25 fluorescent, or the like, are employed to detect events or processes, loaded 

microparticles are spread on a planar substrate, e.g. a glass slide, for examination with 
a scanning system, such as described in International patent applications 
PCT/US91/09217, PCT/NL90/00081, and PCT/US95/01886. The scanning system 
should be able to reproduciblv scan the substrate and to define the positions of each 

30 microparticle in a predetermined region by way of a coordinate system. In 
polynucleotide sequencing applications, it is important that the positional 
identification of microparticles be repeatable in successive scan steps. 

Such scanning systems may be constructed from commercially available 
components, e.g. x-y translation table controlled by a digital computer used with a 

35 detection system comprising one or more photomultiplier tubes, or alternatively, a 
CCD array, and appropriate optics, e.g. for exciting, collecting, and sorting 
fluorescent signals. In some embodiments a confocal optical system may be 
desirable. An exemplary scanning system suitable for use in four-color sequencing is 
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illustrated diagrammatically in Figure 5. Substrate 300, e.g. a microscope slide with 
fixed microparticles, is placed on x-y translation table 302, which is connected to and 
controlled by an appropriately programmed digital computer 304 which may be any of 
a variety of commercially available personal computers, e.g. 486-based machines or 
5 PowerPC model 71 00 or 8100 available form Apple Computer (Cupertino, CA). 

Computer software for table translation and data collection functions can be provided 
by commercially available laboratory software, such as Lab Windows, available from 
National Instruments. 

Substrate 300 and table 302 are operationally associated with microscope 306 

1 0 having one or more objective lenses 308 which are capable of collecting and 

delivering light to microparticles fixed to substrate 300. Excitation beam 310 from 
light source 312, which is preferably a laser, is directed to beam splitter 314, e.g. a 
dichroic mirror, which re-directs the beam through microscope 306 and objective lens 
308 which, in turn, focuses the beam onto substrate 300. Lens 308 collects 

1 5 fluorescence 3 1 6 emitted from the microparticles and directs it through beam splitter 
314 to signal distribution optics 318 which, in turn, directs fluorescence to one or 
more suitable opto-electronic devices for converting some fluorescence characteristic, 
e.g. intensity, lifetime, or the like, to an electrical signal. Signal distribution optics 
318 may comprise a variety of components standard in the art, such as bandpass 
20 filters, fiber optics, rotating mirrors, fixed position mirrors and lenses, diffraction 

gratings, and the like. As illustrated in Figure 5, signal distribution optics 318 directs 
fluorescence 316 to four separate photomultiplier tubes, 330, 332, 334, and 336, 
whose output is then directed to pre-amps and photon counters 350, 352, 354, and 
356. The output of the photon counters is collected by computer 304, where it can be 
25 stored, analyzed, and viewed on video 360. Alternatively, signal distribution optics 
3 1 8 could be a diffraction grating which directs fluorescent signal 318 onto a CCD 
array. 

The stability and reproducibility of the positional localization in scanning will 

determin e, to a large extent, the resolution for separating closely spaced 

30 microparticles. Preferably, the scanning systems should be capable of resolving 

closely spaced microparticles, e.g. separated by a particle diameter or less. Thus, for 
most applications, e.g. using CPG microparticles, the scanning system should at least 
have the capability of resolving objects on the order of 10-100 m. Even higher 
resolution may be desirable in some embodiments, but with increase resolution, the 
35 time required to fully scan a substrate will increase; thus, in some embodiments a 
compromise may have to be made between speed and resolution. Increases in 
scanning time can be achieved by a system which only scans positions where 
microparticles are known to be located, e.g from an initial full scan. Preferably, 
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microparticle size and scanning system resolution are selected to permit resolution of 
fluorescently labeled microparticles randomly disposed on a plane at a density 
between about ten thousand to one hundred thousand microparticles per cm^. 

In sequencing applications, loaded microparticles can be fixed to the surface 
5 of a substrate in variety of ways. The fixation should bfe strong enough to allow the 
microparticles to undergo successive cycles of reagent exposure and washing without 
significant loss. When the substrate is glass, its surface may be derivatized with an 
alkylamino linker using commercially available reagents, e.g. Pierce Chemical, which 
in turn may be cross-linked to avidin, again using conventional chemistries,; to form 

10 an avidinated surface. Biotin moieties can be introduced to the loaded microparticles 
in a number of ways. For example, a fraction, e.g. 10-15 percent, of the cloning 
vectors used to attach tags to polynucleotides are engineered to contain a unique 
restriction site (providing sticky ends on digestion) immediately adjacent to the 
polynucleotide insert at an end of the polynucleotide opposite of the tag. The site is 

15 excised with the polynucleotide and tag for loading onto microparticles. After" 
loading, about 10-15 percent of the loaded polynucleotides will possess the unique 
restriction site distal from the microparticle surface. After digestion with the 
associated restriction endonuclease, an appropriate double stranded adaptor 
containing a biotin moiety is ligated to the sticky end. The resulting microparticles 

20 are then spread on the avidinated glass surface where they become fixed via the 
biotin-avidin linkages. 

Alternatively and preferably when sequencing by ligation is employed, in the 
initial ligation step a mixture of probes is applied to the loaded microparticle: a 
fraction of the probes contain a type lis restriction recognition site, as required by the 

25 sequencing method, and a fraction of the probes have no such recognition site, but 
instead contain a biotin moiety at its non-ligaring end. Preferably, the mixture 
comprises about 10-15 percent of the biotinylated probe. 

In still another alternative, when DNA-loaded microparticles are applied to a 
glass substrate, the DNA may nonspecifically adsorb to the glass surface upon several 

30 hours, e.g. 24 hours, incubation to create a bond sufficiently strong to permit repeated 
exposures to reagents and washes without significant loss of microparticles. 
Preferably, such a glass substrate is a flow cell, which may comprise a channel etched 
in a ulass slide. Preferably, such a channel is closed so that fluids may be pumped 
through it and has a depth sufficiently close to the diameter of the microparticles so 

35 that a monolayer of microparticles is trapped within a defined observation region. 
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Parallel Sequencing 
The tagging system of the invention can be used with single base sequencing 
methods to sequence polynucleotides up to several kilobases in length. The tagging 
system permits many thousands of fragments of a target polynucleotide to be sorted 
5 onto one or more solid phase supports and sequenced simultaneously. In accordance 
with a preferred implementation of the method, a portion of each sorted fragment is 
sequenced in a stepwise fashion on each of the many thousands of loaded 
microparticles which are fixed to a common substrate— such as a microscope slide- 
associated with a scanning system or an image analysis system, such as described 

10 above. The size of the portion of the fragments sequenced depends of several factors, 
such as the number of fragments generated and sorted, the length of the target 
polynucleotide, the speed and accuracy of the single base method employed, the 
number of microparticles and/or discrete regions that may be monitored 
simultaneously; and the like. Preferably, from 12-50 bases are identified at each 

1 5 microparticle or region; and more preferably, 1 8-30 bases are identified at each 
microparticle or region. With this information, the sequence of the target 
polynucleotide is determined by collating the 12-50 base fragments via their 
overlapping regions, e.g. as described in U.S. patent 5,002,867. The following 
references provide additional guidance in determining the portion of the fragments 

20 that must be sequenced for successful reconstruction of a target polynucleotide of a 
given length: Lander and Waterman, Genomics, 2: 231-239 (1988); Drmanac et al, 
Genomics, 4: 1 14-128 (1989); Bains, DNA Sequencing and Mapping, 4: 143-150 
(1993); Bains, Genomics, 11: 294-301 (1991); Drmanac et al, J. Biomolecular 
Structure and Dyriarnics, 8: 1085-1 102 (1991); and Pevzner, J. Biomolecular 

25 Structure and Dynamics, 7: 63-73 (1989). Preferably, the length of the target 

polynucleotide is between 1 kilobase and 50 kilobases. More preferably, the length is 
between 10 kilobases and 40 kilobases. Lander and Waterman (cited above) provide 
guidance concerning the relationship among the number of fragments that are 
sequenced (i.e. the sample size), the amount of sequence information obta in ed from 

30 each fragment, and the probability that the target polynucleotide can be reconstructed 
from the partial sequences without gaps, or "islands." For the present invention, 
maximal polynucleotide sizes that can be obtained for given sample sizes and sizes of 
fragment sequences are shown below: 



35 



-34- 



WO 96/41011 

Size of Sample 



PCT/US96/09513 
Approx. maximai target polynucleotide length 



30 bases/fragment 50 bases/fragment 



1,000 3 kilobases 4 kilobases 

10,000 22 kilobases 32 kilobases 

20,000 40 kilobases 65 kilobases 

30,000 60 kilobases 85 kilobases 

100,000 180 kilobases 300 kilobases 



Fragments may be generated from a target polynucleotide in a variety of ways, 
including so-called "directed" approaches where one attempts to generate sets of 

5 fragments covering the target polynucleotide with minimal overlap, and so-called 
"shotgun" approaches where randomly overlapping fragments are generated. 
Preferably, "shotgun" approaches to fragment generation are employed because of 
their simplicity and inherent redundancy. For example, randomly overlapping 
fragments that cover a target polynucleotide are generated in the following 

10 conventional "shotgun" sequencing protocol, e.g. as disclosed in Sambrook et al 

(cited above). As used herein, "cover" in this context means that every portion of the 
target polynucleotide sequence is represented in each size range, e.g. all fragments 
between 100 and 200 basepairs in length, of the generated fragments. Briefly, starting 
with a target polynucleotide as an insert in an appropriate cloning vector, e.g. phage, 

15 the vector is expanded, purified'and digested with the appropriate restriction enzymes 
to yield about 10-15 \x% of purified insert. Typically, the protocol results in about 
500-1000 subclones per microgram of starting DNA. The insert is separated from the 
vector fragments by preparative gel electrophoresis, removed from the gel by 
conventional methods, and resuspended in a standard buffer, such as TE (Tris- 

20 EDTA). The restriction enzymes selected to excise the insert from the vector 

preferably-leave-compatible-stiGky-ends on-to 

ligated in preparation for generating randomly overlapping fragments. As explained 
in Sambrook et al (cited above), the circularized DNA yields a better random 
distribution of fragments than linear DNA in the fragmentation methods employed 

25 below. After self-ligating the insert, e.g. with T4 ligase using conventional protocols, 
the purified ligated insert is fragmented by a standard protocol, e.g. sonication or 
DNase 1 digestion in the presence of Nto** After fragmentation the ends of the 
fragments are repair, e.g. as described in Sambrook et al (cited above), and the 
repaired fragments are separated by size using gel electrophoresis. Fragments in the 

30 300-500 basepair range are selected and eluted from the gel by conventional means, 
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and ligated into a tag-carrying vector as described above to form a library of tag- 
fragment conjugates. 

As described above, a sample containing several thousand tag-fragment 
conjugates are taken from the library and expanded, after which the tag-fragment 
5 inserts are excised from the vector and prepared for specific hybridization to the tag 
complements on microparticles, as described above. Depending of the size of the 
target polynucleotide, multiple samples may be taken from the tag-fragment library 
and separately expanded, loaded onto microparticles and sequenced. As discussed 
above, the number of doubles selected will depend on the fraction of the tag repertoire 

1 0 represented in a sample. (The probability of obtaining triples—three different 

polynucleotides with the same tag- or above can safely be ignored). As mentioned 
above, the probability of doubles in a sample can be estimated from the Poisson 
distribution p(double)==m 2 e- m /2, where m is the fraction of the tag repertoire in the 
sample. Table VI below lists probabilities of obtaining doubles in a sample for given 

1 5 tag size, sample size, and repertoire diversity. 

Table VIII 



Number of 
words in tag 
from 8 word set 


Size of tag 
repertoire 


Size of 
sample 


Fraction of 
repertoire 
sampled 


Probability of 
double 


7 


2.1 x 10° 


3000 


1.43 x 10-3 


10"° 


8 


1.68 xlO 7 


3 xlO 4 


1.78 x 10" 3 


1.6 xlO" 6 






3000 


1.78 xlO" 4 


1.6 xlO" 8 


9 


1.34x108 


3 xlO 5 


2.24 x 10" 3 


2.5 x 10- 6 






3x 10 4 


2.24 x 10* 4 


2.5 x 10- 8 


10 


1.07 xlO 9 


3x10° 


2.8 x 10- 3 


3.9 x 10- 6 






3 xlO 5 


2.8 x 10" 4 


3.9 x 10" 8 



20 In any case, the loaded microparticles are then dispersed and fixed onto a glass 

microscope slide, preferably via an avidin-biotin coupling. Preferably, at least 15-20 
nucleotides of each of the random fragments are simultaneously sequenced with a 
single base method. The sequence of the target polynucleotide is then reconstructed 
by collating the partial sequences of the random fragments by way of their 

25 overlapping portions, using algorithms similar to those used for assembling contigs, 
or as developed for sequencing by hybridization, disclosed in the above references. 
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Kits for Implen. anting the Method of the Invention 
The invention includes kits for carrying out the various embodiments of the 
invention. Preferably, kits of the invention include a repertoire of tag complements 
attached to a solid phase support. Additionally, kits of the invention may include the 
5 corresponding repertoire of tags, e.g. as primers for amplifying the polynucleotides to 
be sorted or as elements of cloning vectors which can also be used to amplify the 
polynucleotides to be sorted. Preferably, the repertoire of tag complements are 
attached to microparticles. Kits may also contain appropriate buffers for enzymatic 
processing, detection chemistries, e.g. fluorescent or chemiluminescent tags, and the 
1 0 like, instructions for use, processing enzymes, such as ligases, polymerases, 

transferases, and so on. In an important embodiment for sequencing, kits may also 
include substrates, such as a avidinated microscope slides, for fixing loaded 
microparticles for processing. 

15 Identification of Novel Polynucleotides 

in cDNA Libraries 

Novel polynucleotides in a cDN A library can be identified by constructing a 
library of cDNA molecules attached to microparticles, as described above. A large 
fraction of the library, or even the entire library, can then be partially sequenced in 
20 parallel. After isolation of mRNA, and perhaps normalization of the population as 
taught by Soares et al, Proc. Natl. Acad. ScL, 91 : 9228-9232 (1994), or like 
references, the following primer may by hybridized to the polyA tails for first strand 
synthesis with a reverse transcriptase using conventional protocols: 

25 5 ' -mRNA- [A] n -3' 

[T] x9~ [primer site] -GG [W, W, W, C] 9 ACCAGCTGATC-5 ■ 

where [W,W,W,C]9 represents a tag as described above, "ACCAGCTGATC" is an 
optional sequence forming a restriction site in double stranded form, and "primer site" 

"30 is~a~seqreriee^ia^ a primer 

binding site for amplifying polynucleotides of interest by PCR. 

After reverse transcription and second strand synthesis by conventional 
ie;--uiiques, the double stranded fragments are inserted into a cloning vector as 
described above and amplified. The amplified library is then sampled and the sample 

35 amplified. The cloning vectors from the amplified sample are isolated, and the tagged 
cDNA fragments excised and purified. After rendering the tag single stranded with a 
polymerase as described above, the fragments are methylated and sorted onto 
microparticles in accordance with the invention. Preferably, as described above, the 
cloning vector is constructed so that the tagged cDNAs can be excised with an 
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endonuclease, such as Fok I, that will allow immediate sequencing by the preferred 
single base method after sorting and ligation to microparticles. 

Stepwise sequencing is then carried out simultaneously on the whole library, 
or one or more large fractions of the library, in accordance with the invention until a 
5 sufficient number of nucleotides are identified on each cDN A for unique 

representation in the genome of the organism from which the library is derived. For 
example, if the library is derived from mammalian mRNA then a randomly selected 
sequence 14-15 nucleotides long is expected to have unique representation among the 
2-3 thousand megabases of the typical mammalian genome. Of course identification 

10 of far fewer nucleotides would be sufficient for unique representation in a library 
derived from bacteria, or other lower organisms. Preferably, at least 20-30 
nucleotides are identified to ensure unique representation and to permit construction 
of a suitable primer as described below. The tabulated sequences may then be 
compared to known sequences to identify unique cDNAs. 

1 5 Unique cDNAs are then isolated by conventional techniques, e.g. constructing 

a probe from the PCR amplicon produced with primers directed to the prime site and 
the portion of the cDNA whose sequence was determined. The probe may then be 
used to identify the cDNA in a library using a conventional screening protocol. 

The above method for identifying new cDNAs may also be used to fingerprint 

20 mRNA populations, either in isolated measurements or in the context of a 
dynamically changing population. Partial sequence information is obtained 
simultaneously from a large sample, e.g. ten to a hundred thousand, or more, of 
cDNAs attached to separate microparticles as described in the above method. The 
frequency distribution of partial sequences can identify mRNA populations from 

25 different cell or tissue types, as well as from diseased tissues, such as cancers. Such 
mRNA fingerprints are useful in monitoring and diagnosing disease states, e.g. 
International application PCT/US95/21944, which describes the use of express 
sequence tags (ESTs) for the same purpose. 



30 Cycle Sequencing on Microparticles Loaded 

with Sorted Polynucleotides 
Parallel sequencing may also be accomplished in accordance with the 
invention with conventional sequencing techniques that require the generation and 
separation of labeled DNA fragments. In particular, isolated microparticles loaded 
35 with a uniform population of templates may be used to generate labeled extension 
products by cycle sequencing. Cycle sequencing is a well-know variant of the basic 
Sanger approach to DNA sequencing describe fully in the following references: 
Craxton, Methods, Vol. 2 (February, 1991); Wozny, European patent publication 0 



-38- 



WO 96/41011 



PCT/US96/09513 



409 078 A2 (23 January 1991); Fuller, International application PCT/US92/07303; 
and Fuller, International application PCT/US94/03264. Briefly, in a standard 
sequencing reaction mixture, a thermal stable polymerase is employed so that 
repeated extension reactions may be carried out on the same template. This permits 

5 small amounts of template to generate sufficient amounts of extension product for 
detection after separation by electrophoresis. Typically, cycle sequencing comprises 
the steps of (a) providing a sequencing reaction mixture with a template, a primer, 
nucleoside triphosphates, chain-terminating nucleoside triphosphates, and a thermal 
stable DNA polymerase; (b) denaturing the template, (c) annealing the primer to the 

1 0 denatured template, (d) extending the primer to form extension products, and (e) 
repeating steps (b)-(d) until sufficient quantities of extension products are 
accumulated so that they may be detected upon separation. The number of times the 
cycle is repeated depends on many factors, including the amount and quality of 
starting template, the detection system employed, the separation system employed, 

15 and the like. As conventionally practiced, the extension cycle is typically repeated 
from 10 to 100 times; the template amount ranges from as little as a few tens of 
femtomole to several tens of picomole; the denaturation step is carried out by heating 
the reaction mixture to a temperature in the range of 92-95°C; the annealing step 
takes place at a temperature in the range of 35-75°C; and the extension step takes 

20 place at a temperature in the range of 65-85°C with a thermal stable DNA 

polymerase, such as Taq or Vent (available from Perkin-Elmer Corp., Norwalk, CT, 
and New England Biolabs, respectively). 

Tag complements may be prepared on magnetic microparticles as described by 
Albretsen et al, Anal. Biochem., 189: 40-50 (1990), which allows loadings of several 

25 femtomoles of tag complements onto 4.5 jim diameter magnetic beads. Tag 

complements may be attached to the microparticles either by their 5* or 3* ends. If 
attached by 5* ends, then the templates may be sorted via specific hybridization of tags 
at their 3' ends. In this embodiment, the template has a primer complement at its 5* 

end,_as_sho_wn_below: . — — 

30 

3'-[oligonucleotide tagHtemplate] -[primer complement]^' 

The tag complement is then extended the length of the template so that a complement 
of the template is obtained which is covalently attached to the microparticle. The 
35 template is removed by heating and the microparticles are washed. After 

microparticles are separated, e.g. by flow sorting, repeated cycles of annealing 
primers, extension, and denaturation are carried out. 
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If tag complements are attached to the microparticles by their 3 f ends, which 
allows for convenient synthesis directly on the microparticles, the order of the 
oligonucleotide tag and primer complement are reversed, as shown below: 

5 5-loligonucleotide tag]-[template] -[primer complement]^* 

Also, the 5' end of the tag complement is phosphorylated, e.g. using commercially 
available reagents. After specific hybridization via the oligonucleotide tag, a primer 
is annealed to the primer complement at the 3' end of the template and extended with 

1 0 a DNA polymerase lacking 3'-»5' exonuclease activity. The nick left by this 

extension reaction is then ligated and the original template removed by heating. After 
separating microparticles, the cycle sequencing can be carried out as above. 

Separation of loaded microparticles may be carried out by flow sorting, 
wherein suspended microparticles are entrained to pass single file through a nozzle 

1 5 and in a liquid jet which is broken up into a regular series of charged droplets which 
are directed to predetennined target vessels, wells, or other reaction locations on a 
substrate. Microparticles are conveniently detected in the jet by light scatter and the 
magnitude of the scatter is used to determine whether a droplet contains no, one, or 
multiple microparticles. A particularly useful apparatus for such flow sorting and 

20 delivery of sequencing reagent is disclosed in Brennan, International application 
PCT/US94/05896. Once the individual loaded microparticles are distributed to a 
plurality of reaction sites or wells with the appropriate sequencing reagents, the 
collection of reactions can be thermally cycled together to generate extension 
products. After cycling is completed, the extension products are separated by 

25 electrophoresis. Preferably, electrophoretic separation is carried out by capillary 

electrophoresis in a gel-free separation medium, which allows convenient loading and 
rapid separation of the extension fragments. Also, apparatus is available which 
permits detection by four-color fluorescence of a large number of samples 

substantially at the same time, e.g. the type disclosed by Mathies and Huang, Nature, 

30 359:167-169 (1992); Huang et al, Anal. Chem., 64: 2149-2154 (1992); Huang et al, 
Anal. Chem., 64: 967-972 (1992); or the like. Preferably, several thousand cycle 
sequencing reactions are carried at the same time. More preferably, mixtures of 
templates are sorted onto a population of microparticles having a repertoire of 
oligonucleotide tags of between 1000 and 10,000 different types. 

35 

Sorting Multi-locus Probes for Genotvpic Analysis 
Many disease conditions and/or disease susceptibilities are associated with 
complex genetic traits and/or patterns of mutation, e.g. HLA type, mutation pattern of 
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the p53 gene in many cancers, the cystic fibrosis gene, Lesch-Nyhan syndrome, 
Duchenne muscular dystrophy, and the like. Lander et al, Science, 265: 2037-2048 
(1994); Collins, Science, 256:774-779 (1992); Tsui et al, International patent 
application PCT/CA90/00267; Hedrum et al, Biotechniques, 17: 118-129 (1994); 
5 Santamaria et al, International patent application PCT/US92/01675; Chamberlain et 
al, Nucleic Acids Research, 16: 1 1 141-1 1 1 56 (1988); and the like. One approach to 
constructing convenient assays for such complex genetic traits has been to use so- 
called multiplex PCR or multiplex ligation assays, such as described in Chamberlain 
et al (cited above) or in Grossman et al, International patent application 
1 0 PCT/US93/03229. Usually, such techniques call for the simultaneous amplification 
of multiple genetic sequences in the same reaction mixture followed by specific 
detection of sequences of interest. Oligonucleotide tags of the invention can provide 
a simple and convenient means for identifying genetic sequences that are amplified in 
such assays. In its simplest form, this embodiment of the invention may be 
1 5 implement by attaching oligonucleotide tags to PCR primers used in multiplex PCR. 
One primer of a pair carries an oligonucleotide tag and the other primer of the pair 
carries a capture moiety, such as described above, that permits isolation and then 
release of successfully amplified sequences. After release, the sequences are applied 
to solid phase support having a set of tag complements attached at predefined spatial 
20 addresses. The pattern of specific hybridization of the tags is then detected to identify 
the genotype of a sample. 

In a preferred embodiment, PCR is employed to amplify a genetic sequences 
of interest that contains multiple target sites, i.e. multiple sites where mutations or~ 
disease-related sequences occur. Preferably, only two or very few pairs of primers are 
25 used to amplify the target sequence to avoid the difficulties involved with multiplex 
PCR, such as balancing target lengths, primer annealing temperatures, and the like. 
After amplification, specific genotypes are detected in a manner analogous to that 
described in Grossman et al (cited above) and Grossman et al, U.S. patent 5,514,543, 

__whichjeferences.p^ the selection of PCR and ligation reaction 

30 conditions, ligation probe sizes, and the like. In those references, a target sequence is 
similarly amplified, after which a collection of ligation probes are applied in the 
presence of a DNA ligase. The ligation probes consist of two separate sequences both 
complementary to a target potentially present in a sample being analyzed: one is 
attached to a electrophoretic mobility modifier and the other is attached to a 
35 fluorescent label. If the two probes form perfect duplexes with the target sequence in 
the sample they are ligated so that the mobility modifying moiety is now attached to a 
fluorescent label through the ligated sequences complementary to the target The 
components of the mixture are then separated electrophoretically so that the pattern of 
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fluorescent bands on a gel is indicative of the genotype of the target present in the 
sample. As shown in Figure 3, oligonucleotide tags of the invention may be used in 
place of the electrophoretic mobility modifiers and spatial separation can be achieved 
by sorting ligated sequences to particular locations on a solid phase support. 
5 Returning to Figure 3, target sequence (200) is amplified, preferably by PCR, after 
which a collection of ligation probes (206-216) is applied (204) to a denatured 
amplicon. In this embodiment, ligation probes comprise an oligonucleotide tag (206), 
a first sequence (208) complementary to a target sequence, a second sequence (210) 
complementary to the target sequence and contiguous with the first sequence (such 
1 0 that if both are perfectly complementary to the target sequence they are capable of 
being ligated), a tail (212) carrying a signal generating means (214). Signal 
generating means (214) is preferably a fluorescent label. Preferably, the first and 
second sequences of the ligation probes are ligated by a DNA ligase; thus, the 5' end 
of the abutting sequences (216) must be phosphorylated, e.g. via a phosphorylaung 
1 5 reagent described in Urdea et al, U.S. patent 5,332,845. After application of the 
ligation probes and a ligase, probes forming perfectly matched duplexes with the 
target sequence are covalently joined (218 & 220). The probe-target duplexes are 
then denatured and applied (222) to a solid phase support which has tag complement 
attached at well defined spatial locations for every tag t! through t k . After washing 
20 off nonspecifically bound sequences, the spatial locations corresponding to the tag 
complements of the oligonucleotide tags, tj and tj , which were ligated to fluorescent 
labels are illuminated, as shown in Figure 3 by 226 and 228. The pattern of 
illuminated fluorophors on the solid phase support indicates the genotype of the target 
sequence in the sample. Preferably, in this embodiment of the invention there is a 
25 one-to-one correspondence between a tag and a spatial address on the solid phase 
support. In further preference, this embodiment is employed to simultaneously 
identify at least twenty gene targets; and more preferably, it is employed to 
simultaneously detect at least 50 gene targets. 

Generally, this embfldimefltM Reinventi on may he with the f o41owjngsteps_ 

30 for detecting the presence or absence of a plurality of selected target sequences in a 
target polynucleotide: (1) adding to the target polynucleotide a plurality of ligation 
probes, each ligation probe including a first oligonucleotide and a second 
oligonucleotide which are complementary in sequence to adjacent portions of a 
selected one of the target sequences in the target polynucleotide, the first 
35 oligonucleotide having an oligonucleotide tag attached, each oligonucleotide tag 

being selected from the same minimally cross-hybridizing set and each ligation probe 
having a different oligonucletide tag; (2) hybridizing the ligation probes with the 
target polynucleotide; (3) treating the hybridized first and second oligonucleotides 
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under conditions effective to ligate the first and second oligonucleotides whenever the 
first and second oligonucleotides form perfectly match duplexes with adjacent target 
sequences; (4) separating ligated first and second oligonucleotides from unligated first 
and second oligonucleotides; (5) sorting the ligated first and second oligonucleotides 
5 by specifically hybridizing the oligonucleotide tags with their respective 

complements, the respective complements being attached as uniform populations of 
substantially identical oligonucleotides in spatially discrete regions on the one or 
more solid phase supports; and (6) detecting the presence or absence of the selected 
target sequences by the presence or absence of ligated first and second 
1 0 oligonucleotides on the one or more solid phase supports. 



Example 1 

Sortine Multiple Target Polynucleotides Derived from pUC19 
1 5 A mixture of three target polynucleotide-tag conjugates are obtained as 

follows: First, the following six oligonucleotides are synthesized and combined 
pairwise to form tag 1, tag 2, and tag 3: 

S'-pTCGACCfWi) <W2) (w 3 ) (w 4 ) (w 5 ) (w 6 ) (w 7 ) (w 8 ) (W!>A 
20 GG(**) (**) (**) (**) (**) (**) (**) (**) (**)TTCGAp-5' 

Tagl 



25 5 r -pTCGACC(w 6 ) (w 7 ) (w 8 ) (wx) (w 2 ) (w 6 ) (w 4 ) (w 2 ) (wx)A 

GG ( * * ) (**) (**) (**) (**) (**) (**) (**) (**)TTCGAp-5' 



30 



Tag 2 



5 , -pTCGACC(w 3 ) (w 2 ) (wi) (vii) (w 5 ) <w 8 ) (w 8 ) (w 4 ) (w 4 )A 

GG (** ) (**) (**) (**) (**) (**) (**) (**) (**)TTCGAp-5' 



Tag3 



35 



where "p M indicates a monophosphate, the Wj*s represent the subunits define in Table 
n, and the terms "(**)" represent their respective complements. A pUC19 is digested 
with Sal I and Hind m, the large fragment is purified, and separately ligated with tags 
1, 2, and 3, to form pUC19-l, pUC19-2, and pUC19-3, respectively. The three 
40 recombinants are separately amplified and isolated, after which pUC19-l is digested 
with Hind m and Aat I, pUC19-2 is digested with Hind III and Ssp I, and pUC19-3 is 
digested with Hind m and Xmn I. The small fragments are isolated using 
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conventional protocols to give three double stranded fragments about 250, 375, and 
575 basepairs in length, respectively, and each having a recessed 3' strand adjacent to 
the tag and a blunt or 3* protruding strand at the opposite end. Approximately 12 
nmoles of each fragment are mixed with 5 units T4 DNA polymerase in the 
5 manufacturer's recommended reaction buffer containing 33 M deoxycytosine 

triphosphate. The reaction mixture is allowed to incubate at 37°C for 30 minutes, 
after which the reaction is stopped by placing on ice. The fragments are then purified 
by conventional means. 

CPG microparticles (37-74 mm particle size, 500 angstrom pore size, Pierce 

1 0 Chemical) are derivatized with the linker disclosed by Maskos and Southern, Nucleic 
Acids Research, 20: 1679-1684 (1992). After separating into three aliquots, the 
complements of tags 1, 2, and 3 are synthesized on the microparticles using a 
conventional automated DNA synthesizer, e.g. a model 392 DNA synthesizer 
(Applied Biosystems, Foster City, CA). Approximately 1 mg of each of the 

1 5 differently derivatized microparticles are placed in separate vessels. 

The T4 DNA polymerase-treated fragments excised from pUC19-l, -2, and -3 
are resuspended in 50 L of the manufacturer's recommended buffer for Taq DNA 
ligase (New England Biolabs). The mixture is then equally divided among the three 
vessels containing the 1 mg each of derivatized CPG microparticles. 5 units of Taq 

20 DNA ligase is added to each vessel, after which they are incubated at 55°C for 1 5 

minutes. The reaction is stopped by placing on ice and the microparticles are washed 
several times by repeated centrifugation and resuspension in TE. Finally, the 
microparticles are resuspended in Nde I reaction buffer (New England Biolabs) where 
the attached polynucleotides are digested. After separation from the microparticles 

25 the polynucleotide fragments released by Nde I digestion are fluorescently labeled by 
incubating with Sequenase DNA polymerase and fluorescein labeled thymidine 
triphosphate (Applied Biosystems, Foster City, CA). The fragments are then 
separately analyzed on a nondenaturing polyacrylamide gel using an Applied 
Biosystems model 373 DNA sequencer. 

30 

Example 2 

Parallel Sequencing of SV40 Fragments 
A repertoire of 36-mer tags consisting of nine 4-nucleotide subunits selected 
from Table II is prepared by separately synthesizing tags and tag complements by a 
35 split and mix approach, as described above. The repertoire is synthesized so as to 

permit ligation into a Sma I/Hind IE digested M13mpl9. Thus, as in Example I, one 
set of oligonucleotides begins with the addition of A followed by nine rounds of split 
and mix synthesis wherein the oligonucleotide is extended subunit-wise by 3- 
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phosphoramidite derivatived 4-mers corresponding to the subunits of Table II. The 
synthesis is then completed with the nucleotide-by-nucleotide addition of one half of 
the Sma I recognition site (GGG), two Cs, and a S'-monophosphate, e.g. via the 
Phosphate-ON reagent available from Clontech Laboratories (Palo Alto, CA). The 

5 other set of oligonucleotides begins with the addition of three Cs (portion of the Sma 
I recognition site) and two G's, followed by nine rounds of split and mix synthesis 
wherein the oligonucleotide is extended by 3'-phosphoramidite derivatized 4-meFs 
corresponding to the complements of the subunits of Table II. Synthesis is completed 
by the nucleotide-by-nucleotide addition of the Hind HI recognition site and a 5- 

10 monophosphate. After separation from the synthesis supports the oligonucleotides are 
mixed under conditions that permit formation of the following duplexes: 

S'-pGGGCCtwi) (Wi) (w ± ) (Wi) (Wi) (w ± ) (wi) (w ± ) (w ± )A 

CCCGG(**) (**) (**) (**) (**) (**) (**) (**) (**)TTCGAp-5' 

15 

The mixture of duplexes is then ligated into a Sma I/Hind Ill-digested M13mpl9. A 
repertoire of tag complements are synthesized on CPG microparticles as described 
above. 

Next the following adaptor is prepared which contains a Fok I site and 
20 portions of Eco RI and Sma I sites: 

5 1 -pAATTCGGATGATGCATGCATCGACCC 

GCCTACTACGTACGTAGCTGGGp- 5 * 

25 Eco RI Fok I Sma I 

The adaptor is ligated into the Eco Rl/Sma I digested Ml 3 described above. 

Separately, SV40 DNA is fragmented by sonication following the protocol set 
forth in Sambrook et al (cited above). The resulting fragments are repaired using 

30 standard protocols and separated by size. Fragments in the range of 300-500 

basepairs are selected and ligated into the Sma I digested Ml 3 described above to 
form a library of fragment-tag conjugates, which is then amplified. A sample 
containing several thousand different fragment-tag conjugates is taken from the 
library, further amplified, and the fragment-tag inserts are excised by digesting with 

35 Eco RI and Hind m. The excised fragment-tag conjugates are treated with T4 DNA 
polymerase in the presence of deoxycytidine triphosphate, as described in Example I, 
to expose the oligonucleotide tags for specific hybridization to the CPG 
microparticles. 
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After hybridization and ligation, as described in Example I, the loaded 
microparticles are treated with Fok I to produce a 4-nucleotide protruding strand of a 
predetermined sequence. A 10:1 mixture (probe 1 :probe 2) of the following probes 
are ligated to the polynucleotides on microparticles. 



Probe 1 FAM- ATCG6ATGAC 

TAGCCTACTGAGCT 

Probe 2 biotin- ATCCAATGAC 

TAGGTTACTGAGCT 

FAM represents a fluorescein dye attached to the 5'-hydroxyl of the top strand of 
Probe 1 through an aminophosphate linker available from Applied Biosystems 
(Aminolinker). The biotin may also be attached through an Aminolinker moiety and 
optionally may be further extended via polyethylene oxide linkers, e.g. Jaschke et al 
(cited above). 

The loaded microparticles are then deposited on the surface of an avidinated 
glass slide to which and from which reagents and wash solutions can be delivered and 
removed The avidinated slide with the attached microparticles is examined with a 
scanning fluorescent microscope (e.g. Zeiss Axioskop equipped with aNewport 
Model PM500-C motion controller, a Spectra-Physics Model 2020 argon ion laser 
producing a 488 nm excitation beam, and a 520 nm long-pass emission filter, or Wee 
apparatus). The excitation beam and fluorescent emissions are delivered and 
collected, respectively, through the same objective lens. The excitation beam and 
collected fluorescence are separated by a dichroic mirror which directs the collected 
fluorescence through a series of bandpass filters and to photon-counting devices 
corresponding to the fluorophors being monitored, e.g. comprising Hamamatsu model 
9403-02 photomultipliers, a Stanford Research Systems model SR445 amplifier and 
^^^^M^ueX scaler, and digM^B^g^4gfeba s ed 
The computer generates a two dimensional map of the slide which registers the 



15 



20 



25 



30 



35 



positions of the microparticles. 

After cleavage with Fok 1 to remove the initial probe, the polynucleotides on 
the attached microparticles undergo 20 cycles of probe ligation, washing, detection, 
cleavage, and washing, in accordance with the preferred single base sequencing 
methodology described below. Within each detection step, the scanning system 
records the fluorescent emission corresponding the base identified at each 
microparticle. Reactions and washes below are generally carried out with 
manufacturer's (New England Biolabs') recommended buffers for the enzymes 
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15 



employed, unless otherwise indicated. Standard buffers are also described in 
Sambrook et al (cited above). 

The following four sets of mixed probes are provided for addition to the target 
polynucleotides: 



TAMRA- ATCGGATGACATCAAC 

TAGCCTACTGTAGTTGANNN 



10 FAM- ATCGGATGACATCAAC 

TAGCCTACTGTAGTTGCNNN 



ROX- ATCGGATGACATCAAC 

TAGCCTACTGTAGTTGGNNN 

JOE- ATCGGATGACATCAAC 

TAGCCTACTGTAGTTGTNNN 



where TAMRA, FAM, ROX, and JOE are spectrally resolvable fluorescent labels 
20 attached by way of Aminolinker II (all being available from Applied Biosystems, Inc., 
Foster City, California); the bold faced nucleotides are the recognition site for Fok I 
endonuclease, and "N" represents any one of the four nucleotides, A, C, G, T. 
TAMRA (tetramethylrhodamine), FAM (fluorescein), ROX (rhodamine X), and JOE 
(2\7-<iimethoxy-4 , ,5 , -dichlorofluorescein) and their attachment to oligonucleotides is 
25 also described in Fung et al, U.S. patent 4,855,225. 

The above probes are incubated in approximately 5 molar excess of the target — 
polynucleotide ends as follows: the probes are incubated for 60 minutes at 16°C with 
200 units of T4 DNA ligase and the anchored target polynucleotide in T4 DNA ligase 
buffer; after washing, the target polynucleotide is then incubated with 100 units T4 
30 polynucleotide kinase in the manufacturer's recommended buffer for 30 minutes at 

37°C, washed, and again incubated for 30 minutes at 16°C with 200 units of T4 DNA 
igase and the anchored target polynucleotide in T4 DNA liga se buffer. Washi ng is 



accomplished by successively flowing volumes of wash buffer over the slide, e.g. TE, 
disclosed in Sambrook et al (cited above). After the cycle of ligation- 

35 phosphorylation-ligation and a final washing, the attached microparticles are scanned 
for the presence of fluorescent label, the positions and characteristics of which are 
recorded by the scanning system. The labeled target polynucleotide, i.e. the ligated 
complex, is then incubated with 10 units of Fok I in the manufacturer's recommended 
buffer for 30 minutes at 37°C, followed by washing in TE. As a result the target 

40 polynucleotide is shortened by one nucleotide on each strand and is ready for the next 
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cycle of ligation and cleavage. The process is continued until twenty nucleotides are 
identified. 

Example 3 

Construction of a Tag Library 
An exemplary tag library is constructed as follows to form the chemically 
synthesized 9-word tags of nucleotides A, G, and T defined by the formula: 

3 , -TGGC-[4(A,G,T) 9 ]-CCCCp 

where M [ 4 ((A,G,T)9]" indicates a tag mixture where each tag consists of nine 4-mer 
words of A, G, and T; and "p" indicate a 5' phosphate. This mixture is ligated to the 
following right and left primer binding regions: 

15 5'- AGTGGCTGGGCATCGGACCG 5'- GGGGCCCAGTCAGCGTCGAT 

TCACCGACCCGTAGCCp GGGTCAGTCGCAGCTA 

LEFT RIGHT 

20 The right and left primer binding regions are ligated to the above tag mixture, after 
which the single stranded portion of the ligated structure is filled with DNA 
polymerase then mixed with the right and left primers indicated below and amplified 
to give a tag library. 



25 



Left Primer 
AGTGGCTGGGCATCGGACCG 



30 5 1 - AGTGGCTGGGCATCGGACCG- [ 4 < ( A, G , T ) 9 ] -GGGGCCCAGTCAGCGTCGAT 

TCACCGACCCGTAGCCTGGC- [ 4 < (A, G, T) 9] -C CCCGGG TCAGT CGCAG CTA 



CCCCGGGTCAGTCGCAGCTA- 5 1 



35 Right Primer 

The underlined portion of the left primer binding region indicates a Rsr E recognition 
site. The left-most underlined region of the right primer binding region indicates 
recognition sites for Bsp 1201, Apa I, and Eco O 1091, and a cleavage site for Hga I. 
40 The right-most underlined region of the right primer binding region indicates the 
recognition site for Hga I. Optionally, the right or left primers may be synthesized 
with a biotin attached (using conventional reagents, e.g. available from Clontech 
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Laboratories, Palo Alto, CA) to facilitate purification after amplification and/or 
cleavage. 

Example 4 

5 Construction of a Plasmid Library of Tag-Pol vnucleotide 

Conjugates for cDNA "Signature" Sequencing 
cDNA is produced from an mRNA sample by conventional protocols using 
pGGCCCTi5(A or G or C) as a primer for first strand synthesis anchored at the 
boundary of the poly A region of the mRNAs and Ng(A or T)GATC as the primer for 

1 0 second strand synthesis. That is, both are degenerate primers such that the second 
strand primer is present in two forms and the first strand primer is present in three 
forms. The GATC sequence in the second strand primer corresponds to the 
recognition site of Mbo I; other four base recognition sites could be used as well, such 
as those for Bam HI , Sph I, Eco RI, or the like. The presence of the A and T adjacent 

15 to the restriction site of the second strand primer ensures that a stripping and 

exchange reaction can be used in the next step to generate a five-base 5 1 overhang of 
"GGCCC". The first strand primer is annealed to the mRNA sample and extended 
with reverse transcriptase, after which the RNA strand is degraded by the RNase H 
activity of the reverse transcriptase leaving a single stranded cDNA. The second 

20 strand primer is annealed and extended with a DNA polymerase using conventional 
protocols. After second strand synthesis, the resulting cDNAs are methylated with 
CpG methylase (New England Biolabs, Beverly, MA) using manufacturer's protocols. 
The 3* strands of the cDNAs are then cut back with the above-mentioned stripping 
and exchange reaction using T4 DNA polymerase in the presence of dATP and dTTP, 

25 after which the cDNAs are ligated to the tag library of Example 3 previously cleaved 
with Hga I to give the following construct: 



S'-biotin- | primer binding site 1 I TAG 1 1 Apa t 1 1 cDNA ] | 

RsrIIsite Mbo I site 



Separately, the following cloning vector is constructed, e.g. starting from a 
35 commercially available plasmid, such as a Bluescript phagemid (Stratagene, La Jolla, 
CA). 
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primer binding site 
i 



Ppu MI site 

4, 



(plasmid) -5 1 -J 



AAAAGGAGGAGGCCTTGATAGAGAGGACCT GTTTAAAC- 



10 



5 




GGATCC-TCTTCCTCTTCCTCTTCC-3 1 - (plasmid) 
CCTAGG-AGAAGGAGAAGGAGAAGG- 



-CAAATTTG 



t 



t 



Bam HI site 



15 



Pme I site 



The plasmid is cleaved with Ppu MI and Pme I (to give a Rsr Il-compatible end and a 
flush end so that the insert is oriented) and then methylated with DAM methylase. 
The tag-containing construct is cleaved with Rsr II and then ligated to the open 
20 plasmid, after which the conjugate is cleaved with Mbo 1 and Bam HI to permit 

ligation and closing of the plasmid. The plasmid is then amplified and isolated and 
used in accordance with the invention. 
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APPENDIX la 

Exemplary computer program for generating 
minimally cross hybridizing sets 
(single stranded tag/single stranded tag complement) 

Program minxh 

c 

c 

c 

integer*2 subl (6) ,msetl (1000, 6) ,mset2 (1000, 6) 
dimension nbase(6) 

c 
c 

write (*#*)' ENTER SUBUNIT LENGTH ' 
read(*, 100)nsub 
100 format (il) 

open ( 1 , f ile= 1 sub4 . dat 1 , f orm= ' formatted • , stat us= ' new ' ) 

c 
c 

nset=0 

do 7000 ml=l,3 

do 7000 m2=l, 3 

do 7000 m3=l, 3 

do 7000 m4=l,3 
subl(l)=ml 
subl (2) =m2 
subl ( 3 ) =m3 
subl ( 4 ) =m4 

c 
c 

ndiff=3 

c 
c 

c Generate set of subunits differing from 

c subl by at least ndiff nucleotides. 

c Save in msetl. 

c 

c 

do— 90 0 j = 1 , rvsub 

900 msetl (1, j )=subl ( j ) 

c 

c 

do 1000 kl-1,3 

do 1000 k2=l,3 
do 1000 k3=l,3 

do 1000 k4=l,3 

c 
c 

nbase (l)=kl 
nbase (2)=k2 
nbase (3)=k3 
nbase (4 )=k4 
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1200 

c 

c 

c 
c 
c 
c 
c 
c 
c 
c 



1100 

c 
c 

1000 

c 

c 



1325 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



n-0 

do 1200 j=l,nsub 
if <subl ( j ) . eq. 
subl ( j ) .eq. 
subl ( j ) . eq, 
n=n+l 
endif 
continue 



if (n . ge .ndif f ) then 



1 .and. nbase(j) .ne.l .or. 

2 .and. nbase { j ) . ne . 2 .or. 

3 .and. nbase { j ) . ne . 3 ) then 



If number of mismatches 
is greater than or equal 
to ndiff then record 
subunit in matrix mset 



do 1100 i=l, nsub 

msetl ( j j, i)=nbase (i) 
endif 



continue 



do 1325 j 2=1, nsub 
mset2 (1, j2)=msetl (1, j2) 
mset2 (2, j2)=msetl (2, j2) 



c 
c 
c 
c 
c 
c 
c 

c 
c 

1700 



Compare subunit 2 from 
msetl with each successive 
subunit in msetl, i.e. 3, 
4,5, ... etc. Save those 
with mismatches .ge. ndiff 
in matrix mset2 starting at 
position 2. 

Next transfer contents 
of mset2 into msetl and 
start 

comparisons again this time 
starting with subunit 3. 
Continue until all subunits 
undergo the comparisons. 



npass-0 



continue 

kk=npass+2 

npass=npass+l 
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do 



m=npass+2, j j 



1600 



1625 

1500 

c 

c 

c 

c 

c 

c 

c 

c 



2000 



c 
c 



7009 

7008 
7010 



120 
7000 



1500 
n=0 

do 1600 j=l,nsub 

if (msetl (npass + 1 , j ) . eq . 1 . 
rasetl (npass+1, j ) . eq . 2 , 
msetl (npass+1 , j ) .eq.3. 
n»n+l 
endif 
continue 
if (n.ge.ndif f ) then 
kk-kk+1 

do 1625 i=l,nsub 

mset2 (kk, i) =msetl <m, i) 

endif 
continue 



and . mset 1 (m, j ) . ne . 1 . or . 
and . mset 1 (m, j ) . ne . 2 . or . 
and. msetl (iri, j ) . ne . 3) then 



kk is the number of subunits 
stored in mset2 

Transfer contents of mset2 
into msetl for next pass. 



do 2000 k=l,kk 

do 2000 m=l,nsub 

msetl (k,m)=mset2 (k,m) 
if (kk.lt. jj) then 
jj=kk 
goto 1700 
endif 



„nset=nset+l . . - 
write(l,7009) 

format (/) 
do 7008 k=l,kk 

write (1, 7010) (msetl (k,m) /m=l, nsub) 
format (4il) 
write(*, *) 

write (*, 120) kk, nset 

format (lx, •Subunits in set=' , i5, 2x, 'Set No=',i5) 
continue 
close ( 1) 



c 
c 



c 
c 



end 



*** + * + *****★* + * + + ** + *** + + ★ + *■★**** 
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APPENDIX lb 

Exemplary computer program for generating 
minimally cross hybridizing sets 
(single stranded tag/single stranded tag complement) 



Program tagN 

c 

c 

c Program tagN generates minimally cross-hybridizing 

c sets of subunits given i) N — subunit length, and ii) 

c an initial subunit sequence. tagN assumes that only 

c 3 of the four natural nucleotides are used in the tags. 

c 

c 

character*! subl(20) 

integer+2 mset (10000, 20) , nbase{20) 

c 
c 

write (*,*) 'ENTER SUBUNIT LENGTH 1 

read (*, 100) nsub 
100 format (i2) 

c 

• C 

write (*,*)• ENTER SUBUNIT SEQUENCE ' 

read(*,110) (subl (k) , k=l,nsub) 
110 format (20al) 
c 
c 

ndiff=10 

c 
c 

c Let a=l c=2 g=3 & t=4 

c 

c 

do 800 kk=l,nsub 

if (subl(kk) .eq. 'a') then 

mset (1, kk)=l 

endif 

if (subMkkJ.eq.'c 1 ) then 
mset(l,kk)=2 
endif 

if (subl(kk).eq. 'g') then 

msetj likk)«3 

endif 

if <subl(kk) .eq. 'f) then 
mset (l,kk)«4 
endif 

800 continue 

c 

c 

c Generate set of subunits differing from 

c subl by at least ndiff nucleotides, 

c 
c 



c 
c 



do 1000 kl=l,3 
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c 
c 



do 1000 k2=l,3 
do 1000 k3=l,3 
do 1000 k4-l,3 
do 1000 k5«l,3 
do 1000 k6=l,3 
do 1000 k7=l,3 
do 1000 k8-l,3 
do 1000 k9«l,3 
do 1000 kl0«l,3 

do 1000 kll=l,3 
do 1000 kl2=l,3 
do 1000 kl3«l,3 
do 1000 kl4»l,3 
do 1000 kl5=l,3 
do 1000 kl6=l,3 
do 1000 kl7~l,3 
do 1000 kl8=l,3 
do 1000 kl9=l,3 
do 1000 k20-l,3 



nbase 
nbase 
nbase 
nbase 
nbase 
nbase 
nbase 
nbase 
nbase 
nbase 
nbase 
nbase 
nbase 
nbase 
nbase 
nbase 
nbase 
nbase 
nbase 
nbase 



<l)~kl 
<2)«k2 

(3) =k3 

(4) -k4 

(5) =k5 

(6) «k6 
<7)«k7 
(8)°k8 
<9)-k9 

(10) »kl0 

(11) «kll 

(12) -kl2 

(13) -kl3 

(14) =kl4 

(15) =kl5 

(16) =kl6 

(17) =kl7 

(18) -kl8 

(19) -kl9 
<20)«k20 



c 
c 



do 1250 nn=l,jj 



1200 

c 



n=0 

do 1200 j«l,nsub 

if (mset (nn, j) .eq.l .and. 
mset (nn, j) .eq.2 .and. 

Tn<g<=* * (nn f j ) .eg. 3_. anc 

mset (nn, j) .eq.4 .and. 



nbase (j ) .ne. 1 .or. 
nbase(j) .ne.2 .or. 
nbase(j) .ne.3 .or. 



nbase { j ) .ne. 4 ) then 



n-n+1 
endif 
continue 



1250 

c 

c 



if (n.lt.ndiff) then 

goto 1000 

endif 
continue 



jj=jj+l 

write (*, 130) (nbase (i) , i=l,nsub) , j j 
do 1100 i=l,nsub 
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mset ( j j , i)«nbase (i) 
1100 continue 
c 
c 

1000 continue 

c 

c 

write (*,*) 
130 format (lOx, 20 (lx, il) , 5x, i5) 

write(*,*> 

write(*, 120) jj 
120 format (lx, •Number of words=\i5) 

c 
c 

end 

c ******************************************** 



************************* 



******************* 
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APPENDIX Ic 

Exemplary computer program for generating 
minimally cross hybridizing sets 
(double stranded tag/single stranded tag complement) 

Program 3tagN 

c 

c 

c Program 3tagN generates minimally cross-hybridizing 

c sets of duplex subunits given i) N — subunit length, 

c and ii) an initial homopurine sequence. 

c 

c 

character+1 subl(20) 

integer*2 mset { 10000, 20) , nbase(20) 

c 
c 

write (*,*) 'ENTER SUBUNIT LENGTH • 

read (*, 100) nsub 
100 format (i2) 
c 
c 

write ( * , * ) • ENTER SUBUNIT SEQUENCE a & g only" 

read(*,110) (subl ( k) , k-1, nsub) 
110 format (20al) 

c 
c 

ndiff=10 

c 
c 

c Let a=l and g=2 

c 

c 

do 800 kk=l,nsub .... 

if (subl(kk) .eq. 'a') then 

mset (l,kk)-l 

endif 

if (subl(kk) .eq. 'g' ) then 
mset (1, kk)«2 
endif 

BOO continue 

c 

c 



do 1000 kl«l,3 
do 1000 k2«l,3 
do 1000 k3«l,3 
do 1000 k4=l,3 
do 1000 k5=l,3 
do 1000 k6«l,3 
do 1000 k7=l,3 
do 1000 k8=l,3 
do 1000 k9=l,3 
do 1000 kl0«l,3 

do 1000 kll=l,3 
do 1000 kl2=l,3 
do 1000 k!3«l,3 
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c 
c 



c 
c 



1200 

c 

c 



do 1000 kl4=l,3 
do 1000 kl5=l,3 
do 1000 kl6=l,3 
do 1000 kl7=l,3 
do 1000 kl8=l,3 
do 1000 kl9=l,3 
do 1000 k20=l,3 



nbase (1) 
nbase (2) 
nbase (3) 
nbase (4 ) 
nbase (5) 
nbase (6) — . 
nbase (7)= 
nbase (8) = 
nbase (9) " 
nbase (10) 
nbase (11) 
nbase (12) 
nbase (13) 
nbase (14) 
nbase (15) 
nbase (16) 
nbase (17) 
nbase (18) 
nbase (19) 
nbase (20) 



kl 
k2 
k3 
k4 
k5 
k6 
k7 
k8 
*9 
=kl0 
«kll 
~kl2 
=kl3 
«kl4 
=kl5 
=kl6 
=kl7 
=kl8 
«kl9 
-k20 



do 1250 nn=l, jj 
n=0 

do 1200 j=l,nsub 
if (mset(nn, j) 

1 mset (nn, j) , 

2 mset(nn,j), 

3 mset (nn, j) , 
n«n+l 
endif 
continue 



eq. 
eq. 
eq. 
eq. 



. and . nbase ( j ) . ne . 1 . or . 
.and. nbase ( j ) .ne. 2 .or. 
.and. nbase(j) .ne.3 .or. 
. and . nbase ( j ) . ne . 4 ) then 



1250 
c 



if (n.lt.ndiff ) then 

goto 1000 

endif 
continue 



JJ-jJ+ 1 



1100 
c 

1000 
c 

130 



120 

c 

c 



write ( * , 130 ) { nbase ( i ) , i=l , nsub ) , j j 
do 1100 i~l,nsub 

mset ( jj, i)«nbase(i) 
continue 

continue 
write (*,*) 

format (lOx, 20 (lx, il) , 5x, i5) 
write {*,*) 
write (*, 120) jj 

format (lx, 'Number of words^iS) 



end 
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SEQUENCE "LISTING 



(1) GENERAL INFORMATION: 



(i) APPLICANT: Sydney Brenner 



(ii) TITLE OF INVENTION: Oligonucleotide Tags for Sorting and 
Identification 



(iii) NUMBER OF SEQUENCES: 16 



(iv) CORRESPONDENCE ADDRESS: 

(A) ADDRESSEE: Stephen C. Macevicz, Spectragen, Inc. 

(B) STREET: 3832 Bay Center Place 

(C) CITY: Hayward 

(D) STATE: California 

(E) COUNTRY: USA 

(F) ZIP: 94545 



(v) COMPUTER READABLE FORM: 

(A) MEDIUM TYPE: 3.5 inch diskette 

(B) COMPUTER: IBM compatible 

CO OPERATING SYSTEM: Windows 3.1 
(D) SOFTWARE: Microsoft Word 5.1 



(vi) CURRENT APPLICATION DATA: 

(A) APPLICATION NUMBER: 
<B) FILING DATE: 
<C) CLASSIFICATION: 

(vii) PRIOR APPLICATION DATA: 

(A) APPLICATION NUMBER: PCT/US95/12791 

(B) FILING DATE: 12-OCT-95 



(vii) PRIOR APPLICATION DATA: 

(A) APPLICATION NUMBER: 08/478,238 

(B) FILING DATE: 07-JUN-95 

(vii) PRIOR APPLICATION DATA: 

(A) APPLICATION NUMBER: 08/485,105 

(B) FILING DATE: 07-JUN-95 

(vii) PRIOR APPLICATION DATA: 

(A) APPLICATION NUMBER: P CT/US95/12791 
(B) FILING~DATE1 12 -OCT- 9 5 

(viii) ATTORNEY/AGENT INFORMATION: 

(A) NAME: Stephen C. Macevicz 

(B) REGISTRATION NUMBER: 30,285 

(C) . REFERENCE / DOCKET NUMBER: cbd4wo 



(ix) TELECOMMUNICATION INFORMATION: 

(A) TELEPHONE: (510) 670-9365 

(B) TELEFAX: (510) 670-9302 



(2) INFORMATION FOR SEQ ID NO: 1: 
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(i) SEQUENCE CHARACTERISTICS: 

{A) LENGTH: 38 nucleotides 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS : single 
{ D) TOPOLOGY: linear 

(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 1: 



GAGGATGCCT TTATGGATCC ACTCGAGATC CCAATCCA 



(2) INFORMATION FOR SEQ ID NO: 2: 



(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 26 nucleotides 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 2: 



AATTCGGATG ATGCATGCAT CGACCC 



(2) INFORMATION FOR SEQ ID NO: 3: 



(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 14 nucleotides 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 3: 



TAGCCTACTG AGCT 



(2) INFORMATION FOR SEQ ID NO; 4: 



(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 16 nucleotides 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: double 

(D) TOPOLOGY: linear 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 4: 



ATCGGATGAC ATCAAC 



(2) INFORMATION FOR SEQ ID NO: 5: 
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(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 11 nucleotides 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: double 

(D) TOPOLOGY: linear 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 5: 



ACCAGCTGAT C 11 
(2) INFORMATION FOR SEQ ID NO: 5: 

(1) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 11 nucleotides 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 5: 

CTAGTCGACC A 11 

(2) INFORMATION FOR SEQ ID NO: 6: 

(i) SEQUENCE CHARACTERISTICS : 

<A) LENGTH: 11 nucleotides 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 
<D) TOPOLOGY: linear 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 6: 

NRRGATCYNN N 11 



-{2T~ INFORMATION-FOR-SEQ-1D NO: TT 



(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 22 nucleotides 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 7: 

GGGTCGATGC ATGCATCATC CG 22 
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(2) INFORMATION FOR SEQ ID NO: 8: 



(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 10 nucleotides 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 8: 



ATCGGATGAC 



(2) INFORMATION FOR SEQ ID NO: 9: 



(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 10 nucleotides 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 9: 



ATCNNNNNAC 



(2) INFORMATION FOR SEQ ID NO: 10: 



(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 14 nucleotides 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 



1 xi ) SEQUENCE^DESCRIPT ION :~ SEQ TD~NO: 10 : 



TCGAGTNNNN NGAT 



(2) INFORMATION FOR SEQ ID NO: 11: 



(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 16 nucleotides 

(B) TYPE: nucleic acid 
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(C) STRANDEDNESS : single 

(D) TOPOLOGY : linear 

(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 11: 

ATCGGATGAC ATCAAC 16 



(2) INFORMATION FOR SEQ ID NO: 12: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 20 nucleotides 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 12: 

NNNAGTTGAT GTCATCCGAT 20 



(2) INFORMATION FOR SEQ ID NO: 13: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 20 nucleotides 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 13: 

NNNCGTTGAT GTCATCCGAT 20 



(2) INFORMATION FOR SEQ ID NO: 14: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 20 nucleotides 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 14: 

NNNGGTTGAT GTCATCCGAT 20 
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(2) INFORMATION FOR SEQ ID NO: 15: 



(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 20 nucleotides 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS : single 

(D) TOPOLOGY: linear 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 15: 



NNNTGTTGAT GTCATCCGAT 20 



(2) INFORMATION FOR SEQ ID NO: 16: 



(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 37 nucleotides 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: double 

(D) TOPOLOGY: linear 

(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 16: 



NNNNNGGATG NNNNNNNNNN NNNTNNNNNN NNNNNNN 
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I claim: 



1 . A method of sorting polynucleotides from a population of polynucleotides onto one 
or more solid phase supports, the method comprising the steps of: 
5 (a) attaching an oligonucleotide tag from a repertoire of tags to each polynucleotide 

in a population of polynucleotides such that each oligonucleotide tag from the repertoire is 
selected from the same minimally cross-hybridizing set; 

(b) sampling the population of polynucleotides such that substantially all different 
polynucleotides in the population have different oligonucleotide tags attached; and 
1 0 (c) sorting the polynucleotides from the population by specifically hybridizing the 

oligonucleotide tags with their respective complements, the respective complements being 
attached as uniform populations of substantially identical oligonucleotides in spatially 
discrete regions on the one or more solid phase supports. 

15 2. The method of claim 1 wherein said oligonucleotide tag is single stranded. 

3. The method of claim 2 wherein said complement of said oligonucleotide tag is single 
stranded. 

20 4. The method of claim 3 wherein said oligonucleotide tag consists of a plurality of 
subunits, each subunit consisting of an oligonucleotide of 3 to 9 nucleotides in length and 
each subunit being selected from the same minimally cross-hybridizing set. 

5. The method of claim 4 wherein said repertoire of said oligonucleotide tags contains 
25 at least 1 000 of said oligonucleotide tags. 

6. The method of claim 5 wherein each of said oligonucleotide tags differs from every 
other oligonucleotide tag of said minimally cross-hybridizing set by at least three 

nucleotides 

30 

7. The method of claim 6 wherein each of said oligonucleotide tags has a length in the 
range of from 12 to 60 nucleotides. 



8. The method of claim 7 wherein said one or more solid phase supports are 
35 microparticles. 

9. The method of claim 8 wherein said microparticles have a diameter in the range of 
from 5 to 40 um. 
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10. The method of claim 9 wherein said repertoire contains at least 10,000 of said 
oligonucleotide tags. 

5 11. The method ofclaim 7 wherein each ofsaid one or more solid phase supports is a 
planar substrate having a plurality of spatially discrete surface regions, each of such surface 
regions having attached thereto a uniform population of substantially identical said 
complements. 

10 12. The method ofclaim 1 1 wherein each ofsaid spatially discrete surface regions has 
an area in the range of from 10 to 1000 jim 2 . 

13. The method ofclaim 3 wherein each ofsaid oligonucleotide tags consists of an 
oligonucleotide having a length between 12 and 30 nucleotides. 



15 



14. The method ofclaim 13 wherein said repertoire contains at least 100 ofsaid 
oligonucletide tags. 



15. The method ofclaim 14 wherein each ofsaid oligonucleotide tags differs from every 
20 other oligonucleotide tag by at least six nucleotides. 

16. The method ofclaim 1 5 wherein each ofsaid oligonucleotide tags consists of an 
oligonucleotide having a length between 15 and 24 nucleotides. 

25 17. The method ofclaim 16 wherein said repertoire contains at least 1000 ofsaid 
oligonucletide tags. 

18. The method ofclaim 1 wherein said oligonucleotide tag is double stranded. 



"30 19^ The method ofclaim 18 wherein said oligonucleotide tag has a first strand and a 
second strand such that the first strand is homopyrimidine and the second strand is 
homopurine. 

20. The method of claim 19 wherein said complement is homopyrimidine. 

21 . The method ofclaim 20 wherein said oligonucleotide tag consists of a plurality of 
subunits, each subunit consisting of an oligonucleotide of 4 to 10 base pairs in length and 
each subunit being selected from the same minimally cross-hybridizing set. 
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22. The method of claim 21 wherein said repertoire of said oligonucleotide tags contains 
at least 100 of said oligonucleotide tags. 

5 23. The method of claim 22 wherein each of said oligonucleotide tags differs from every 
other oligonucleotide tag of said minimally cross-hybridizing set by at least three base pairs. 

24. The method of claim 23 wherein each of said oligonucleotide tags has a length in the 
range of from 12 to 60 nucleotides. 

10 

25. The method of claim 24 wherein said one or more solid phase supports are 
microparticles. 

26. The method of claim 25 wherein said microparticles have a diameter in the range of 
15 from 5 to 40 uin. 

27. The method of claim 26 wherein said repertoire contains at least 1000 of said 
oligonucleotide tags. 

20 28. The method of claim 24 wherein each of said one or more solid phase supports is a 
planar substrate having a plurality of spatially discrete surface regions, each of such surface 
regions having attached thereto a uniform population of substantially identical said 
complements. 

25 29. The method of claim 28 wherein each of said spatially discrete surface regions has 
an area in the range of from 10 to 1000 pjn 2 . 

30. The method of claim 20 wherein each of said oligonucleotide tags consists of an 

oligcmudeotideJuLving-aJengt^ 

30 

3 1 . The method of claim 30 wherein said repertoire contains at least 1 0 of said 
oligonucletide tags. 

32. The method of claim 3 1 wherein each of said oligonucleotide tags differs from every 
35 other oligonucleotide tag by at least six base pairs. 

33. The method of claim 32 wherein each of said oligonucleotide tags consists of an 
oligonucleotide having a length between 15 and 24 base pairs. 
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34. The method of claim 33 wherein said repertoire contains at least 100 of said 
oligonucletide tags. 

5 35. A method of identifying a population of mRNA molecules, the method comprising 
the steps of: 

forming a population of cDNA molecules from the population of mRNA molecules 
such that each cDNA molecule has an oligonucleotide tag attached, the oligonucleotide tags 
being selected from the same minimally cross-hybridizing set; 
10 sampling the population of cDNA molecules such that substantially all different 

cDNA molecules have different oligonucleotide tags attached; 

sorting the cDNA molecules by specifically hybridizing the oligonucleotide tags with 
their respective complements, the respective complements being attached as uniform 
populations of substantially identical complements in spatially discrete regions on one or 
1 5 more solid phase supports; 

determining the nucleotide sequence of a portion of each of the sorted cDNA 
molecules; and 

identifying the population of mRNA molecules by the frequency distribution of the 
portions of sequences of the cDNA molecules. 

20 

36. The method of claim 35 wherein said oligonucleotide tag is single stranded. 

37. The method of claim 36 wherein said complement of said oligonucleotide tag is 
single stranded. 

25 

38. The method of claim 37 wherein said oligonucleotide tag consists of a plurality of 
subunits, each subunit consisting of an oligonucleotide of 3 to 9 nucleotides in length and 
each subunit being selected from the same minimally cross-hybridizing set. 



30 39. The method of claim 38 wherein substantially all of said different oligonucleotide 
tags differ from one another by having at least two different said subunits. 

40. The method of claim 39 wherein said one or more solid phase supports are 
microparticles and wherein said step of sorting said cDNA molecules onto the microparticles 
35 produces a subpopulation of loaded microparticles and a subpopulation of unloaded 
microparticles. 
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41 . The method of claim 40 further including a step of separating said loaded 
microparticles from said unloaded microparticles. 

42. The method of claim 41 further including a step of repeating said steps of sampling, 
sorting, and separating until a number of said loaded microparticles is accumulated is at least 
10,000. 

43. The method of claim 42 wherein said number of loaded microparticles is at least 
100,000. 

44. The method of claim 43 wherein said number of loaded microparticles is at least 
500,000. 

45. The method of claim 41 further including a step of repeating said steps of sampling, 
sorting, and separating until a number of said loaded microparticles is accumulated is 
sufficient to estimate the relative abundance of a cDN A molecule present in said population 
at a frequency within the range of from 0.1% to 5% with a 95% confidence limit no larger 
than 0.1% of said population. 

20 46. A method of detecting the presence or absence of a plurality of selected target 
sequences in a target polynucleotide, the method comprising the steps of: 

adding to the target polynucleotide a plurality of ligation probes, each ligation probe 
including a first oligonucleotide and a second oligonucleotide which are complementary in " 
sequence to adjacent portions of a selected one of the target sequences in the target 
25 polynucleotide, the first oligonucleotide having an oligonucleotide tag attached, each 

oligonucleotide tag being selected from the same minimally cross-hybridizing set and each 
ligation probe having a different oligonucletide tag; 

hybridizing the ligation probes with the target polynucleotide; 

treatmg_theJiybridized_fiist_and.see^^ — 

30 to ligate the first and second oligonucleotides whenever the first and second oligonucleotides 
form perfectly match duplexes with adjacent target sequences; 

separating ligated first and second oligonucleotides from unligated first and second 
oligonucleotides; 

sorting the ligated first and second oligonucleotides by specifically hybridizing the 
35 oligonucleotide tags with their respective complements, the respective complements being 
attached as uniform populations of substantially identical oligonucleotides in spatially 
discrete regions on the one or more solid phase supports; and 



10 
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detecting the presence or absence of the selected target sequences by the presence or 
absence of Jigated first and second oligonucleotides on the one or more solid phase supports. 

47. The method of claim 46 wherein said second oligonucleotide includes a capture 
5 moiety and a reporter moiety. 

48. The method of claim 47 wherein said one or more solid phase supports is a planar 
substrate having a plurality of spatially addressable discrete surface regions, each of such 
surface regions having attached thereto a uniform population of substantially identical said 

1 0 complements. 

49. The method of claim 48 wherein said capture moiety is a biotin. 

50. The method of claim 49 wherein each of said spatially addressable discrete surface 
1 5 regions has an area in the range of from 1 0 to 1 000 um 2 . 

5 1 . The method of claim 50 wherein said plurality is at least twenty. 

52. A method of identifying polynucleotides, the method comprising the steps of: 

20 ( a ) Providing a repertoire of oligonucleotide tags having a complexity substantially 

greater than the complexity of a population of polynucleotides, the oligonucleotide tags 
being selected from the same minimally cross-hybridizing set; 

(b) ligating the repertoire of oligonucleotide tags to the population of 
polynucleotides to form a mixture of tag-polynucleotide conjugates; 
25 (c) sampling the mixture of tag-polynucleotide conjugates to obtain a subpopulation 

of tag-porynucleotide conjugates such that substantially all tag-polynucleotide conjugates 
having different polynucleotides have different oligonucleotide tags; and 

(d) identifying the polynucleotides of the subpopulation of tag-polynucleotide 
conjugates by specifically hybridizing the oligonucleotide tags with their respective 



30 complements. 

53. The method of claim 52 wherein said oligonucleotide tag is single stranded. 

54. The method of claim 53 wherein said complement of said oligonucleotide tag is 
35 single stranded. 
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55. The method of claim 54 wherein said oligonucleotide tag consists of a plurality of 
subunits, each subunit consisting of an oligonucleotide of 3 to 9 nucleotides in length and 
each subunit being selected from the same minimally cross-hybridizing set. 

5 56. The method of claim 55 wherein said complexity of said repertoire is at least ten 
times said complexity of said population of polynucleotides. 

57. The method of claim 56 wherein each of said oligonucleotide tags has a length in the 
range of from 12 to 60 nucleotides. 

10 

58. The method of claim 57 wherein said one or more solid phase supports are 
microparticles. 

59. The method of claim 57 wherein each of said one or more solid phase supports is a 
1 5 planar substrate having a plurality of spatially discrete surface regions, eachof such surface 

regions having attached thereto a uniform population of substantially identical said 
complements. 

60. The method of claim 54 wherein each of said oligonucleotide tags has a length in the 
20 range of from 12 to 30 nucleotides. 

61 . The method of claim 60 wherein each of said oligonucleotide tags differs from every 
other oligonucleotide tag by at least six nucleotides. 

25 62. The method of claim 61 wherein each of said oligonucleotide tags consists of an 
oligonucleotide having a length between 15 and 24 nucleotides. 

63. The method of claim 62 wherein said complements of said oligonucleotide tags are 

attached as uniform populations of substantially identical oligonucleotides in spatjally 

30 addressable regions on one or more solid phase supports. 

64. The method of claim 52 wherein said oligonucleotide tag is double stranded. 

65. The method of claim 64 wherein said oligonucleotide tag has a first strand and a 
35 second strand such that the first strand is homopyrimidine and the second strand is 

homopurine. 

66. The method of claim 65 wherein said complement is homopyrimidine. 
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67. The method of claim 66 wherein said oligonucleotide tag consists of a plurality of 
subunits, each subunit consisting of an oligonucleotide of 4 to 10 base pairs in length and 
each subunit being selected from the same minimally cross-hybridizing set. 

5 

68. The method of claim 67 wherein each of said oligonucleotide tags consists of an 
oligonucleotide having a length between 12 and 30 base pairs. 

69. A repertoire of oligonucleotide tags, the repertoire being selected from the group 
1 0 consisting of oligonucleotides of the form: 

SiS 2 S 3 ... S n 

wherein each of S] through S n are subunits consisting of an oligonucleotide having a length 
from three to nine nucleotides and being selected from a minimally cross-hybridizing set; 
15 n is in the range of from 4 to 1 0, with the proviso that the oligonucleotide tags of the 

repertoire have a length in the range of from 12 to 60 nucleotides or base pairs; and 
wherein there are at least a hundred oligonucleotide tags in the repertoire. 

70. A repertoire of oligonucleotide tags having a length in the range between 1 2 and 30 
20 nucleotides or base pairs, the oligonucleotide tags of the repertoire being selected from the 

same minimally cross-hybridizing set and the oligonucleotide tags of the repertoire differing 
from one another by at least six nucleotides or base pairs. 

71. The repertoire of claim 70 having a complexity of at least one hundred. 

72. The repertoire of claim 71 having a complexity of at least one thousand. 
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